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EXCKALOE  INTERACTIONS  Ili 
EARTH  PETALS. 

R.  J.  Elliott 

Clarendon  Laboratory, 

Oxford,  England. 

ABSTRACT 

The  complex  magnetic  orderings  observed  in  the  rare 

( 1 ) 

earth  metals  may  be  understood  phenomenologically'-  '  in 

terms  of  an  effective  exchange  interaction  J(R)  between  the 

<  * 

localised  4f  electrons  a  distance  R  apart  which  is  of  long 
range  so  that  the  Fourier  transform 

J(a)  =  21  J(s)e12.-S 

has  a  maximum  for  £  at  some  general  point  in  the  Brillouin 
Zone.  The  anisotropy  of  the  systems  can  be  understood  in 
terms  of  the  crystalline  electric  field  which  influences  the 
orbital  motion  of  the  f  electrons.  In  the  second  half  of 
the  rare  earth  series  the  crystal  field  effects  are  com¬ 
parable  to  the  exchange  energy,  in  the  first  half  they  are 
somewhat  larger.  The  detailed  form  of  J(cl)  and  of  the 
anisotropy  terms  cam  be  studied  through  the  spin-wave 

spectra  which  are  ver,,  complex  in  these  complicated  magnetic 
(2) 

structures' . '  The  details  of  the  magnetic  ordering  are 
found  by  neutron  diffraction  to  change  with  temperature  T 
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and  with  aiioying^^ partly  because  of  tne  anisotropy,  t.i  t 
also  because  the  effective  J(_q)  appears  to  change  with  f. 

The  origin  of  the  exchange  interaction  is  believed  to 
lie  in  indirect  exchange  through  the  conduction  electrons^ 
A  nearly  free  electron  calculation  gives  the  correct  * 

/  e  \ 

essential  properties  of  J(£)^'  although  the  actual  Fermi 
surface  must  be  very  different  to  this^^  and  the  effect  is 
enhanced  by  the  interelectronic  interaction^^.  When 
magnetic  ordering  sets  in  the  conduction  electrons  see  an 
exchange  field  of  a  new  symmetry  which  causes  gaps  to 
appear  in  the  one-electron  energy  spectrum.  These  gaps 
follow  new  zone  boundaries  which  tend  to  form  near  the 
Fermi  surface.  The  subsequent  distortion' of  the  Fermi 
surfhee  gives  rise  to  strong  anomalies  in  electrical  trans¬ 
port  properties  like  the  resistivity,  and  in  optical 
absorption*.  It  also  causes  a  change  in  J(p.)  due  to  the 
redistribution  of  electrons. Furthermore  the  scattering 
of  the  conduction  electrons  off  the  spin  disorder  causes  a 
modification  of  the  conduction  electron  polarisation  which 
changes  These  mechanisms  give  a  reasonable 

explanation  of  the  observed  variations. 

Thus  the  combination  of  indirect  exchange  and  crystal 
field  effects  appears  to  account  qualitatively  for  the 
complex  properties  of  the  rare  earth  metals.  A  more 
detailed  test  of  the  theory  awaits  more  experimental  data 


» 
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on  J(,g.)  and  on  the  precise  form  of  the  Fermi  surface. 
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SO ME  THEORETICAL  ASPECTS 
OF  HEAVY  RARE  EARTH  METALS 

Takeo  Nagamiya 

Faculty  of  Engineering  Science 
Osaka  University 
Toyonaka,  Japan 


The  band  structure  for  the  conduction  electrons  and  its 
consequence  on  the  magnetic  and  electric  properties  of  heavy 
rare  earth  metals  are  discussed.  .  A  brief  description  of 
the  theory  of  the  magnetization  processes  of  a  helical  spin 
arrangement,  in  particular  structure  changes  from  helix  to 
fan  and  then  to  ferromagnetic  alignment  with  increasing  mag¬ 
netic  field,  is  presented,  the  temperature  being  assumed  to 
be  finite  and  a  sixfold  anisotropy  within  the  easy  plane  be¬ 
ing  included.  Further,  spin  waves  in  a  helical  arrangement 

i  * 

and  in  a  fan  arrangement  are  discussed;  uniform  modes  and 
their  frequencies  are  particularly  considered. 

1  -  Introduction 

The  theoretical  understanding  of  the  magnetic  and  elec¬ 
tric  properties  of  rare-earth  metals  is  at  present  made 
through  two  lines  of  approach.  (a)  One  assumes  for  the 
conduction  electrons  the  picture  of  the  free  electrons  per¬ 
turbed  by  the  periodic  lattice  potential,  or  of  electrons  in 
oand3  that  would  follow  with  the  single  orthogonalized  plane 
wave  (30PW)  functions.  Then,  one  considers  the  exchange 
and  Coulomb  Interactions  between  the  conduction  electrons 
and  the  Af  shells.  On  this  basis,  one  can  understand 
roughly,  but  consistently,  the  ferromagnetism  of  Gd,  the 
helical  ordering  of  the  Af  moments  in  heavier  metals,  the 
ceriod  of  this  ordering,  its  dependence  on  atomic  number  and 

PRECEDING  PAGE  BLANK 
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temperature,  the  electrical  resistivity  as  a  function  of 
temperature  and  ordering,  and  so  on  (Kasuya1  1964).  How¬ 
ever.  reliable  calculations  of  the  band  structure  and  ex¬ 
perimental  determinations  of  the  Permi  surface  are  not  yet 
available,  which  puts  a  difficulty  in  the  way  of  a  further 
theoretical  study.  (b)  One  assumes  certain  constants  of 
the  indirect  exchange  coupling  among  the  4f  moments  and  also 
certain  forms  of  the  anisotropy  energy  for  each  4f  moment 
which  depend  on  the  4f  shell  electronic  configuration.  One 
can  then  discuss  the  various  moment  orderings  observed  or 
predictable  for  rare-earth  metals.  One  can  also  discuss 
the  temperature  and  magnetic-field  dependence  of  the  order- 

2  "5 

ing,  as  well  as  the  spin  wave  spectra  (Kaplan,  Elliott, 

4  cj 

Miwa  and  Yosida  1961;  Nagamiya  et  al.  1962,  Kitano  and 
Nagamiya^  1964;  Cooper  et  al.^’®  1962,  63).,  However,  the 
values  of  the  theoretical  parameters  have  to  be  sought  ulti¬ 
mately  in  the  interaction  between  the  4f  shells  and  the  con¬ 
duction  electrons,  and  for  this  purpose  one  has  again  to 
know  the  band  structure.  Also,  experimental  data  are  not 
sufficiently  available  to  determine  these  parameters. 

In  this  note,  the  theory  of  the  type  of  (a),  which  is 
inevitably  qualitative,  will  be  briefly  discussed,  and  then 
the  writer's  theory  of  the  type  of  (b)  of  the  magnetization 
processes  in  heavy  rare-earth  metals  will  be  described, 
though  not  fully.  Besides,  the  modes  ana  frequencies  of 
the  spin  waves  in  helical  and  fan  spin-structures,  particu¬ 
larly  of  those  which  are  to  be  resonant  to  external,  homo¬ 
geneous,  oscillating  magnetic  field,  will  be  presented,  in  a 
more  complete  and  correct  form  than  those  published  earlier 
by  Cooper  et  al. 

2.  Conduction  Band  and  Magnetic  and  Electric  Properties 

Q 

In  the  early  theory  of  Yosida  and  tfatabe^  the  conduc¬ 
tion  electrons  were  assumed  to  be  free  at  the  starting  point 
and  then  the  f-s  exchange  interaction  was  introduced  as  a 
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perturbation.  When  the  f-shell  moments  are  aligned  heli¬ 
cally,  there  appear  energy-gap  planes  in  the  k-space  for  the 
conduction  electrons.  If  the  wave  vector,  Q,  of  this  heli¬ 
cal  order  is  such  that  a  certain  set  of  the  gap  planes  come 
to  tontact  with  the  Fermi  surface,  then  the  helical  order  is 
stabilized.  It  was  found  that  such  planes  are  those  which 
perpendicularly  bisect  the  vectors  K^-j+Q  drawn  from  the 
origin.  Here  the  index  111  refers  to  the  orthorhombic  unit 
cell  for  the  h.c.p.  lattice,  and  there  are  altogether  six 
equivalent  indices  of  this  kind.  Q  was  found  to  be  in  the 
c-direction  and  to  correspond  to  a  period  of  approximately 
seven  hexagonal  layers.  Thus,  their  theory  appeared  to 
work  nicely. 

If  we  are  allowed  to  treat  both  the  lattice  potential 
and  the  exchange  potential  as  perturbations,  then  we  shall 
find  that  the  energies  arising  from  these  two  are  additive, 
as  far  as  the  second  order  perturbation  energies  are  con¬ 
cerned;  thus,  the  theory  of  Yosida  and  Watabe  will  be  valid¬ 
ated.  In  the  theory  of  SOPW,  there  is  a  theorem  that  the 
kinetic  energy  increase  of  an  electron  in  a  SOPW  state  due 
to  the  core  functions,  added  to  a  plane  wave  function  for 
the  purpose  of  orthogonalization,  approximately  cancel  the 
potential  energy  decrease  due  to  the  core  potential.  If 
this  cancellation  is  nearly  perfect,  then  we  have  an  energy 
spectrum  of  a  nearly  free  electron,  and  we  can  validate  the 
Yosida-Watabe  theory. 

However,  observations  of  the  electronic  specific  heats 
of  some  rare-earth  metals,  as  well  as  of  Y  and  Sc,  show  that 
they  have  values  seven  times  as  large  as  that  one  would  ex¬ 
pect  with  the  assumption  of  free  electrons.  The  conduction 
electrons  must,  therefore,  be  strongly  perturbed  by  the 
lattice  periodic  potential.  Fig.1  shows  the  ( iTo)  cross- 
section  of  the  Fermi  sphere  for  free  electrons  and  of  the 
Brillouin  zone  boundaries.  An  expected  distortion  of  the 
Fermi  surface  due  to  the  lattice  periodic  potential  is  also 
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shown.  In  the  reduced  zone  scheme,  this  distorted  Fermi 
surface  will  be  as  shown  in  Fig. 2.  One  will  see  by  a  con¬ 
sideration  in  three  dimensions  that  there  are,  in  addition, 
electron-filled  hemispheres  on  the  top  and  bottom  surfaces 
of  the  hexagonal  cylinder  (Fig. 3)  and  holes  at  the  midpoints 
of  the  vertical  edges  of  this  cylinder  (not  shown  in  the 
figure ) . 

We  assume  that  the  cancellation,  mentioned  before,  is 
less  perfect,  thus  the  distortion  of  the  Fermi  surface 
greater,  for  heavier  metals;  in  particular,  the  portion  of 
the  Fermi  surface  at  the  central  part  of  Fig. 2  is  assumed  to 
be  flatter  for  heavier  metals.  If  one  considers  this  por¬ 
tion  only  and  studies  the'  q-dependent  susceptibility (q  being 
the  wave  vector  of  a  hypothetical 


Fig. 1  Pig. 3 
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q-dependent  susceptibility  has  a  maximum  at  q=Q,  where  Q  is 
a  vertical  vector  shown  in  Pig. 2.  This  would  mean  that  a 
helical  order  having  a  wave  vector  Q  stabilizes  itself  by 
producing  gap  planes  which  contact  the  Fermi  surface  as 
shown  by  horizontal  broken  lines  also  in  Pig. 2  (one  plane  is 
produced  by  one  kind  of  spin,  the  other  by  the  other  kind  of 
spin).  The  heavier  the  metal,  the  flatter  will  be  this 
portion  of  the  Permi  surface,  hence  greater  will  be  the 
value  of  Q.  This  is  in  accordance  with  the  observation  of 
the  wave  number  of  the  ordered  arrangement  in  heavy  rare- 
earth  metals  immediately  below  the  N6el  temperature.  With 
decreasing  temperature,  the  magnitude  of  the  energy  gap  will 
increase,  and  the  flat  part  of  the  Permi  surface  contacted 
by  the  gap  plane  will  be  replaced,  more  or  less  completely, 
by  this  plane.  Then,  in  order  to  minimize  the  total  ener¬ 
gy,  the  increase  in  the  energy  gap  will  cause  Q  diminish, 
the  opposite  Permi  surface  being  pushed  away;  this  may  be 
seen  by  a  simple  calculation  for  a  one -dimensional  system 
having  a  parabolic  energy  spectrum.  Thus,  Q  decreases  with 
decreasing  temperature,  as  boserved. 

When  one  has  a  sinusoidal  order,  as  in  Sr  and  Tm,  both 
the  upper  and  lower  gap  planes  appear  for  each  kind  of  spin, 
and  they  will  behave  as  boundaries  for  the  electron-filled 
region.  The  value  of  Q,  therefore,  will  not  change  with 
decreasing  temperature.  Of  course,  the  actual  situation 
may  not  be  so  simple.  In  particular,  the  appearance  of 
harmonics,  in  the  periodic  order  will  certainly  play  an  im¬ 
portant  role  in  fixing  the  period  to  an  integral  number  of 
'layers'0  (seven  layers  in  the  cases  of  Sr  and  Tm). 

We  have  so  far  neglected  other  portions  of  the  Permi  sur¬ 
face.  Including  the  small  electron  surfaces  shown  in  Pig. 3 
and  the  hole  surfaces,  they  will  contribute  a  decreasing 
function  of  |qi  to  the  q-dependent  susceptibility  when  q 
varies  in  the  c-direction.  In  Gd,  in  which  the  maximum  at 
q=Q  in  the  aforethought  susceptibility  may  be  assumed  to  be 
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the  least,  this  contribution  might  smear  it  out  or  make  it 
the  second  maximum.  If  this  is  the  case,  we  expect  ferro¬ 
magnetism  for  Gd. 

The  anisotropy  of  the  Permi  surface  of  the  sort  consi¬ 
dered  above  is  further  supported  by  the  fact  that,  in  the 
paramagnetic  temperature  range,  the  resistivity  parallel  to 
the  c-axis  is  smaller  than  that  perpendicular  to  the  c-axis 
(less  average  effective  mass  in  the  c-direction  than  perpen¬ 
dicular  to  it).  This  anisotropy  in  the  resistivity  in¬ 
creases  with  increasing  atomic  number,  as  expected.  A  hump 
in  the  resistivity  parallel  to  c  versus  T  curve,  observed 
near  Tc  in  Gd  and  near  T^  in  some  heavier  metals,  may  be  in¬ 
terpreted  as  being  due  to  short  range  spin  order  that  corres¬ 
ponds  to  Q.  Further  anomalies  at  lower  temperatures  can 
naturally  be  interpreted  as  being  due  to  long  range  order 
and  the  consequent  decrease  in  the  free  Fermi  surface.11’12 

Kasuya1  discusses  quantitatively  all  these  features  and 
others,  but  the  present  writer  feels  that  we  cannot  go  much 
further  at  present  than  the  qualitative  discussion  given  in 
this  section. 

3*  Structure  Changes  of  a  Helical  Spin  Arrangement  by  Ex¬ 
ternal  Magnetic  Field 

It  has  been  known  since  some  years  ago  that  a  helical 
spin  ordering,  when  subjected  to  a  magnetic  field  applied 
in  the  easy  plane  in  which  the  spin  vectors  rotate,  changes 
discontinuously  to  a  fan-like  ordering  when  the  field  ex¬ 
ceeds  a  certain  critical  value  and  further,  above  another 

critical  field  the  latter  transforms  to  a  parallel  alignment. 

1  3 

The  theory  x*as  put  forward  by  Herein  and  Keriel  and  Nasa- 

miya  et  al.  They  confined  themselves  to  T=0.  The  theory 

for  finite  T  and  for  various  modifications  of  the  helical 
ordering  was  then  worked  out  by  Kitano  and  Nagamiya. 

We  shall  discuss  here  only  the  case  of  helix-fan-f erro 
transitions.  The  temperature  may  be  finite  and  an  aniso- 


11 


troc.y  energy  of  sixfold  symmetry  in  tne  easy  plane  may  be 
included.  The  applied  field  be  in  the  easy  direction  in 
the  nlane .  The  mathematics  Dresented  here  is  different 
from  that  of  reference  6. 

Denoting  by  H*  (vector)  the  effective  field  acting  on 
the  nth  spin,  by  Sdni  the  thermal  average  of  the  spin  compo¬ 
nent  Sni  (i=x,y,z).  and  by  wQ  the  anisotropy  Hamiltonian  of 
the  nth  spin  which  is  a  function  of  Sn,  namely  wn=w( snh  we 
can  write  the  molecular  field  equation  as 

Stfni=  TrCS^xptkT)-1  (H*«S-wn)  ]./Tr  [exp(kT ) (H^-S~wn)3^ 

x  (1) 

We  want  to  solve  this  equation. 

To  study  the  stability  of  the  ferromagnetic  alignment 
along  the  applied  field,  whose  direction  will  be  taken  as  x, 

we  put  <3‘nx=:cro+crnx  and  assume  tha‘t  ^nx’  ^ny?  and  anz  to  de 
small  (the  easy  plane  is  taken  as  the  xy-plane).  Here  6^ 

is  assumed  to  satisfy  the  equations 


Sdn=  Tr[S  exp(kT)”1  (H^S  -wri)]/Tr(exp(kT)_l  (H«S  -w  )  ],  (2) 


-1  ,TT*r 


For  the  y 


'“(Tx  ^crx  n/J’  v'~' 

HS=Zm  2JmnSa0  +  H  =  2J(0)Sfl0  +  H- 

J _  is  the  exchange  constant  between  and  S„ .  For  the  y 

mn  °  m  n  ' 

and  z 'components ,  we  have 

2JmnS5ny'  HL=ZB  2JmnSanz’  <4) 

For  the  x-component  we  put 

Hnx=  Ho  +  Hnx’  then  K  ^nS°nx‘  ™ 

We  want  to  linearize  £q .  ( 1  )  with  respect  to  0.'  ,  <3  ,  a 

#  *  ”*  ny  nz 

and  thus  to  Hnx>  Hny,  Snz. 

Before  doing  that,  we  specialize  the  anisotropy  Hamil¬ 
tonian  as  follows : 

Hn=  »Snz2  -  <G^6S4)t(Snx+iSny)6+(Snx-iSny)6'].  (6: 

In  the  classical  limit,  we  may  put  Snx=  Scoscp^ ,  S  =  Ssinyn< 

when  the  spin  is  in  the  xy-plane;  then,  the  second  term  of 
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(6)  reduces  to  SS~<f^ 


Also,  when  expanded  in  powers  of 


S~Snx  and  Sny,  this  term  can  be  written  to  the  lowest  order 

const  +  -p(S-Snx)  +  gGSny  . 

Now,  linearizing  JSq.(l),  we  may  obtain  equations  of  the 

form  *  *  * 

Stf  =  H  B,  Sff  =  HO,  S<5'  =  H'_A*  (7) 

ny  ny  9  nz  nz  ’  nx  nx 

where  A,  B,  C  are  functions  of  HQ/kT  divided  by  kT.  They 
represent  certain  susceptibilities  along  the  x,  y,  z  direc¬ 
tions,  respectively.  It  can  be  shown,  as  might  be  expected, 
that  for  D  greater  than  G  (both  positive)  B  is  the  largest, 

C  medium,  and  A  the  smallest.  This  means  that  the  ferro¬ 
magnetic  alignment  in  the  x-direction  is  the  least  stable 


for  spin  deviations  in  the  y-dlrection.  Therefore,  we  put 
in  the  first  equation  of  (7)  0yexp(iq*Rn) ,  where  Rn  is 

the  position  of  the  nth  spin.  Then  we  find,  using  the 
first  equation  of  (4), 


<y=  2J(q)BGy,  'where  J(q )=  J^expf iq-B^) .  .  (8) 

Here  R  =  R  -  R  . 
mn  m  n 

Prom  (8)  one  may  see  that  the  ferromagnetic  alignment 
in  the  x-direction  becomes  unstable  when  the  applied  field 
becomes  less  than  the  field,  H  .  at  which  2J(q)B=1.  Below 

C 

this  field  an  oscillation  In  the  y-component  of  the  spins 
will  set  in  with  a  wave  vector  q.  The  quantity  B  can  be 
calculated  to  be  SO'0/(Hq+  2ScTqG),  which  is  a  decreasing 
function  of  H.  Thus,  the  actual  critical  field,  below 
which  the  ferromagnetic  alignment  becomes  unstable,  namely, 
the  largest  possible  field  consistent  with  2J(q)B=1,  must 
correspond  to  the  maximum  of  J(q).  The  value  of  q  which 
gives  this  maximum  is  exactly  the  wave  vector,  Q,  of  the 
helical  ordering.  (We  assume  that  the  function  J(q)  is  not 
affected  by  changes  in  the  spin  ordering,  which  may  not  be 
true  in  the  actual  case.)  Thus,  putting  q=Q,  we  obtain 

Hc=  2S<J0CJ(Q)-J(°)-G]. 


(9) 
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-ere  we  used  Sa  .  ( 3)  to  rewrite  <2J(Q)3=1.  (syig  was  assumed 
to  be  unity  in  the  oresent  notation;  to  convert  to  the  ordi¬ 
nary  notation,  K  is  replaced  by  gp^H.  ) 

We  have  seen  that  the  spin  vectors  oscillate  in  the  xy~ 
plane  with  the  wave  vector  Q  in  the  .field  below  H  .  This 
oscillation  is  a  static  one,  namely  an  oscillation  in  space. 
The  spin  vectors  form  a  fan.  Its  amplitude,  hence  also  the 
susceptibility  of  the  fan  structure  for  variations  of  the 
static  field,  H,  in  the  x-direction  can  be  calculated.  For 
T=0,  the  susceptibility  is  (g2p.g2  times) 

1 /[ 3 J ( Q ) -2 J ( 0 ) - J{ 2Q ) -35G ] •  (10) 

It  may  be  noted  that  this  quantity  can  become  negative  when 
G  increases.  If  this  happens,  the  above  argument  fails. 

One  expects  in  this  case  a  transition  of  the  first  kind 
(discontinuous  transition)  between  the  fan  and  ferromagnetic 
structures  at  a  field  greater  than  that  given  by  (9)*  In 
fact,  a  discontinuous  transition  was  found  by  Kitano  and 
Nagamiya  in  a  certain  range  of  G  with  complicated  calcula¬ 
tions,  including  numerical  computation.  One  has  even  a 
discontinuous  transition  between  the  helical  and  ferromag¬ 
netic  structures  without  an  intermediate  fan  structure  when 
G  is  sufficiently  large,  while  for  smaller  G  double  transi¬ 
tions  helix-fan-ferro  are  predicted.  Koehler  observed  two 
intermediate  fan  structures  in  Ho  in  a  certain  range  of 
temperature,  but  this  is  not  yet  understood  theoretically. 

4.  Spin  Waves  in  the  Helical  and  Fan  Structures 

Spin  waves  in  the  helical  spin  arrangement,  studied 
1  4 

first  by  Yoshimori,  have  the  following  characteristic. 

The  mode  for  q=0  is  a  mere  rotation  of  the  system  about  the 
axis  of  the  helix,  provided  that  there  is  no  anisotropy  in 
the  plane  of  spin  rotation.  Its  frequency  is  therefore  0. 
For  long  waves,  the  spin  wave  frequency  is  proportional  to 
the  wave  number.  It  then  drops  to  a  minimum  near  q=Q  when 
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the  out-of -plane  anisotropy  constant  is  small.  When  the 
latter  is  zero,  then  the  minimum  is  also  zero.  This  can  be 
seen  by  considering  the  mode  for  q=Q.  This  mode  is  such 
that  the  spin  at  R  and  the  spin  at  R„  oscillate  with  a 
phase  difference  of  Q*Rflin»  so  that  the  spins  pointing  oppos¬ 
ite  to  each  other  oscillate  with  a  phase  difference  of  %; 

the  plane  spanned  by  the  spin  vectors  performs  a  tilting 

p 

motion.  Assuming  an  anisotropy  energy  DS  ,  its  frequency 

HZ 

can  be  calculated  to  be 

2S{D[J(Q)-J(0)-J(2Q)]}'/2. 

There  are  actually  two  degenerate  modes  corresponding  to  q=Q 
and  q=-Q,  which  rotate  opposite  to  each  other.  If  a  field 
is  applied  in  the  plane,  in  the  x-direction,  the  frequency 
will  split  into  two,  one  corresponding  to  the  cos-mode  in 
which  the  spins  in  the  x-direction  have  the  largest  ampli¬ 
tude  and  the  other  to  the  sin-mode  in  which  the  spins  per¬ 
pendicular  to  the  x-direction  have  the  largest  amplitude. 

The  frequency  of  the  q=0  mode  has  been  shown  to  increase  as 
2 

H  .  When  the  field  is  above  Hc,  so  that,  we  have  a  ferro¬ 
magnetic  alignment,  the  mode  for  q-0  is  Just  the  ordinary 
resonance  mode  having  a  finite  frequency.  A  peculiar  fact 
is  that  the  frequency  of  the  ferromagnetic  spin  wave  with  a 
wave  number  Q  vanishes  at  H  ,  which  could  he  imagined  from 

V 

that  a  fan  structure  with  Q  starts  to  develop  at  this  field. 

In  the  fan  region  for  H<H  ,  bith  the  modes  for  q=0  and 
q-Q  have  interesting  features.  Q  and  -Q  are  mixed  to  give 
a  cos-mode  and  a  sin-mode,  of  which  the  former  has  a  vanish¬ 
ing  frequency.  In  the  cos-mode,  the  spin  vectors  at  the 
center  of  the  fan  oscillate  in  the  plane  with  the  largest 
amplitude  and  the  spins  at  the  edges  of  the  fan  are  at  rest 
(Pig. 4a).  When  one  sees  the  vectors  opposite  to  the  field 
direction,  one  sees  a  sinusoidal  curve  for  the  locus  of  the 
top  of  the  vectors;  with  the  oscillation  of  the  vectors  this 
sinusoidal  curve  will  shift  up  and  down  as  illustrated  in 


■v 


^**Jvrr 
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Fi.-r.Ab.  Thus,  evidently,  the  frequency  must  vanish. 

In  the  sin-mode,  the  frequency  vaiies  as  (H  -  H)1//2, 
namely,  it  is  proportional  to  the  angular  amplitude  of  the 
fan.  Denoting  the  latter  by  <*,  we  can  write  the  frequency 

^Q,sin=  SQ4i^J(Q)'-2J(0)-J(2Q)-35G+Nx][D-(5/6)G]j  1/2, 

provided  that  «i  is  small  and  the  transition  between  the  fan* 

and  the  ferromagnetic  alignment  is  of  the  second  order.  Here 

N  is  the  demagnetizing  coefficient  in  the  x-direction, 

X  o 

multiplied  by  (gpB)  N,  N  being  the  number  of  atoms  in  unit 

volume.  This  mode  is  illustrated  in  Fig. 5a  and  5b. 

The  mode  for  q=0  is  such  that  the  spin  vectors  describe 

equal  ellipses  in  space,  so  that  the  fan  oscillates  as  a 

rigid  body  in  the  easy  plane  and,  with  a  phase  difference  of 

ti/2,  oscillates  in  the  direction  perpendicular  to  the  easy 

plane  with  some  bending,  like  an  ordinary'  fan  oscillates 

when  people  use  it  to  produce  breeze.  Its  frequency  is 

given  by  /p 

••  ±*Q=  2S[(A0+  B0)(A0-  B0)V/d , 

where 

Aq+  Bq=  J( Q ) -J (0 ) -  |-<*2CJ(Q)-J(0)+35G]-(2S)"1  (Hc«  H) 

+ 

Aq-  Bq=  J(Q)-J(0)  +  (D+  gG)  -  G  +  jf<*2[J(Q)-.J(0)-5G] 

-  (2Sr’(H0-  H)  +  i-Nz. 
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Experimental  Study  of  the  Excitations  in  Rare-Earth  Metals 

A.  R.  MACKINTOSH 

Institute  for  Atomic  Research  and  Department  of  Physics 
Iowa  State  University,  Ames,  Iowa 


Those  properties  of  the  rare  earth  metals  which  are  of  particular 
interest  to  the  solid  state  physicist  are  primarily  determined  by  the 
presence  of  a  well-defined  magnetic  moment  on  the  ion.  The  ordering  of 
these  localized  moments  below  characteristic  temperatures  can  have  a 

profound  effect  on  the  conduction  electrons,  while  the  spin  wave  exci- 

( 

tat  ions  from  the  magnetic  ground  state  are  of  great  importance  in 
determining  the  thermal,  transport  and  magnetic  properties  of  the  metal. 
Considerable  progress  has  been  made  recently,  both  theoretical.ly  and 
experimentally,  in  understanding  these  excitations  and  this  paper  is  an 
attempt  to  review  recent  experimental  developments. 

A  direct  experimental  study  of  the  Fermi  surface  in  the  rare  earth 
metals  is  rendered  difficult  by  the  technical  problems  of  preparing  pure 
single  crystals  and  the  more  fundamental  difficulty  that  the  application 
of  a  magnetic  field,  upon  which  the  most  successful  techniques  depend, 
may  modify  the  magnetic  structure  and  hence  the  Fermi  surface  itself. 

The  study  of  the  angular  correlation  of  the  photons  emitted  when  positrons 
annihilate  with  the  electrons  is  a  relatively  crude  technique,  but  it  does 
not  suffer  from  these  limitations  and  recent  measurements  on  single  crystals 
of  holmium^  have  shown  clearly  that  the  Fermi  surface  is  highly  anisotropic. 
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2 

in  agreement  with  the  anomalously  high  electronic  specific  heat  and 
recently  band  structure  calculations^.  It  is  likely  that  such  experi¬ 
ments,  together  with  calculations,  will  soon  lead  to  a  reasonably 
accurate  picture  of  the  Fermi  surface.  In  particular,  it  should 
readily  be  possible  to  test  the  widely  accepted  hypothesis  that  the 
Fermi  surfaces  of  all  the  hexagonal  close  packed  rare  earths  are  closely 
similar.  The  modification  of  the  electronic  distribution  by  the  magnetic 

ordering  may  be  observed  by  this  technique  and  also  by  infrared  absorption 

.  4 

measurements  . 

Less  direct  information  on  the  conduction  electrons  and  their  inter¬ 
action  with  the  spin  waves  is  given  by  the  study  of  transport  properties. 
The  temperature'*  and  magnetic  field^  dependence  of  the  electrical  resis¬ 
tivity  support  the  view  that  the  Fermi  surface  is  anisotropic,  and 
measurements  along  the  c  axis  give  an  indication  of  the  considerable 
modification  of  the  electronic  distribution  by  the  superzone  planes.^ 

These  planes  apparently  affect  the  resistivity  in  the  basal  plane  very 
little,  so  that  changes  which  occur  at  magnetic  transitions  reflect  pri¬ 
marily  the  difference  in  electron-magnon  scattering  and  give  some  infor¬ 
mation  about  the  spin  waves.  The  electron-magnon  interaction  apparently 

O 

makes  a  large  contribution  to  the  low  temperature  thermoelectric  power, 
through  the  mechanism  of  magnon  drag.  Some  characteristics  of  the  spin 

wave  spectra  can  also  be  determined  from  the  temperature  dependence  of 

,  g  10 

the  magnetization  and  heat  capacity 


Nt*,**-V'*' 


— 
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The  most  direct  and  powerful  method  of  measuring  the  magnon  dispersion 
curves  is  by  inelastic  neutron  scattering,  but  again  experimental  diffi¬ 
culties  appear  in  the  form  of  the  high  thermal  neutron  capture  cross  section 
of  most  of  the  rare  earths.  Experiments  are  currently  being  performed  on 
inelastic  neutron  scattering  in  terbium^,  however,  and  it  seems  that  it 
should  also  be  possible  to  make  a  detailed  study  of  holmium.  Magnetic 

1 2 

resonance  also  gives  directly  the  spin  wave  energies  at  some  wavevectors 
and  this  technique  should  provide  a  valuable  adjunct  to  neutron  scattering, 
particularly  when  the  latter  is  technically  difficult. 

It  appears  that  a  good  description  of  the  Fermi  surface  and  spin  wave 
spectra  of  some  of  the  rare  earths  will  soon  be  available,  and  it  will  then 
be  possible  to  make  quantitative  calculations  of  those  physical  properties 
which  depend  upon  the  electrons,  the  magnons  and  their  interaction. 
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THE  INDIRECT  EXCHANGE  INTERACTION 
IN  THE  RARE  EARTH  METALS* 

Frederick  Specht 

William  Marsh  Rice  University 
Houston,  Texas 


An  exact  calculation  is  given  for  the  indirect  exchange 
interaction  in  the  rare  earth  metals.  With  a  free  electron 
gas  model  for  the  conduction  band  the  isotropic  and  first 
order  non-isotroplc  terms  are  obtained,  and  it  Is  shown 
that  the  non-isotropic  tern  gives  sizeable  contributions  to 
the  interaction  between  two  ions.  An  estimate  is  made  of 
the  contribution  of  the  anisotropy  to  the  crystalline  ex¬ 
change  energy  for  the  ferromagnetic  phase.  The  two  ferro¬ 
magnetic  ordering  patterns  with  moments  along  the  c-axls  and 
in  the  hexagonal  plane  are  considered,  and  an  energy  differ¬ 
ence  in  agreement  with  experiment  is  found  to  be  of  the 
order  of  ten  per  cent. 

In  the  rare  earth  metals  the  principal  interaction 

between  the  4f  shell  electrons  on  different  ions  is  the  in- 

« 

direct  exchange  interaction  via  the  conduction  band.  The 
methods  used  thus  far  in  calculating  this  interaction  have 
involved  some  type  of  approximation  on  the  functional 
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dependence  of  J(k>k')j  the  exchange  Integral  between  the 
iff  shell  electrons  and  a  conduction  electron.  The  approx- 
filiation  used  most  frequently  is  to  replace  J(k,k*)  by  a 
constant.  With  a  spherically  symmetric  conduction  band 
the  leading  term  of  the  interaction  (isotropic)  then  has 
the  radial  dependence  given  by  the  well  known  Rudderman- 
Kittel  ^-Kasuya  ^  ^-Yosida  ^  (RKKY)  function  f  (x  cos(x) 

-  sin(x) )  J  where  x  =  2kpR .  (kp  is  the  Fermi  wave  vec¬ 
tor  of  the  conduction  band  and  R  is  the  distance  between 
the  ions.)  A  somewhat  better  approximation  is  to  set 

J(k,k')  =  J(k  -  k’),  Using  this  approximation,  Kaplan  and 

m 

Lyons'  '  calculated  the  isotropic  and  first  order  non-iso¬ 
tropic  terms  of  the  interaction  and  showed  that  the  radial 
dependence  of  both  terms  differs  from  the  (RKKY)  function. 

We  have  evaluated  the  interaction  exactly  and  obtain  results 
quite  different  from  all  previous  calculations.  Here  we 
give  an  outline  of  the  calculation  and  briefly  consider  the 
effect  of  the  resulting  interaction  on  the  magnetic  order¬ 
ing. 


According  to  the  usual  expression  from  second  order 
perturbation  theory  the  effective  exchange  Hamiltonian  be- 
tween  ions  a  and  b  takes  the  form' 

Hab  =  *  %  (t) 

x  Jaa'  tu: a,  jh ) 


.**-£■  v. 


where  Jaa  *  denotes  the  integral  operator 

J  F d3A  L  exp [i  (2- A')-  (A  -A  ) 

t.  Oi  “* 

and  where  J(k,k’,  X  )  is  defined  by  the  relation 

</■*. '/  J r(A,4'jtJI/*'>  » 

x  m  n(n )  ylm^A2£i^ll 

lx,  ~  £./ 

'  -A. 

R(r)  is  the  4f  shell  radial  function,  and  (aai,£al)  de¬ 
note  the  (spin, orbital)  angular  momentum  operators  of  the 
i’th  electron  on  ion  a.  denotes  a  sum  over  all  the 

4f  shell  electrons  of  the  appropriate  ion.  It  1s  conven- 

*p 

lent  to  define  the  z  axis  in  the  direction  of  Bab- 

The  isotropic  and  first  order  non-isotropic  terms  are 
obtained  by  using  the  spherical  expansions  of  /x-xI/*'/ 
and  exp(ik«r)  and  retaining  only  the  terms  of  the  expansion 

of  exp(ik*r)  with  zero  magnetic  quantum  number  to  second 

a 

order.  J(k,k’,  £  )  then  becomes  a  linear  combination  of 
the  Irreducible  tensor  operators  ^/00(j£  )  and 

•Jb 

Substituting  this  expression  for  J(k,k,,</)  into  Eq,(l) 
and  neglecting  the  term  containing  ^20</al>  ^20 (/  bj^ 
one  obtains  the  exchange  Hamiltonian 

Hal,  =  Z.  sa  ■  S b  Q 


1 
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where  Q  is  of  the  form 

Q  =  JAA'  fa  [J^ (A,  A1)  J/o;  (A', A) 

X  (cntA$,  cnJig  ) [ £fita  Mat )  ^  Jfvo 

()  =  0,1) 

and  where  iAi;  are  integrals  of  the  form 

[Jr,  lJf,  txJfx  [K  (r, )  R  (n  )  Jt  (Ar,) 

*  Jm(A'fz)  (KeVrS+’J] 

*r  / 

The  evaluation  of  the  integrals  7ll'  J./iV  J 
Is  straight  forward.  The  integrations  over  ky  ,k\^  , 

kQ  and  k'g  are  trivial.  The  integrations  over  k  and 
k1  can  be  performed  by  expanding  the  four  Bessel  functions 
in  a  Taylor  series.  Each  term  of  the  expansion  leads  to 
integrals  of  the  form 

XOQ  OO 

X  iMtx  (x-a.ir'Ax  f  XKa*.K(x-a.V'Ax 

V  Jo 

with  m  odd  and  n  even.  With  a  convergence  factor  these 

integrals  become  [{  T  fe)  a01-1  cos  (a  >7  and£(  -  7T/2)  an”^ 

S  ^  (a  \J  respectively.  After  regrouping  the  terms  of  the 
expansion,  one  obtains  _ 

/”  ~x3~ 

j£[  +af  mFj  f  -J  zXv{AP)iMuF) 

X  -  2  sA  F  R  4 1 
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Using  this  procedure,  we  have  evaluated  the  isotropic 
and  first  order  non-isotropic  terras.  We  evaluated  the 
functions  (kF)  and  their  derivatives  by  numerical 

Integration  with  the  Hartree-Fock  4f  shell  functions  given 
by  Freeman  and  Watson^)  #  The  resulting  Hamiltonian  is 
given  by 

*  [  Zfzo  (u )  +  y*0(U) 7/  it) 

where  ^2q(L)  =  (L-n)2  -  L(L+l)/3.  is  the  unit  vector 
between  the  ions,  and  oC  is  defined  by  the  operator  equiv¬ 
alence  ,  .  „  ,  . 

f  $ ,  Wfc  <4>-* 

within  the  manifold  of  the  4f  shell  angular  functions. 

S  and  L  are  the  total  spin  and  orbital  angular  momentum 
operators  of  the  4f  shell.  S,  L  and  are  given  in 
Table  I  for  the  last  half  rare  earths. 

From  the  values  of  a  and  L  it  is  clear  that  the  non¬ 
isotropic  term  gives  a  sizeable  contribution  compared  with 
the  isotropic  term.  Furthermore  both  the  isotropic  and 
the  non-isotropic  functions  differ  considerably  from  the 
(RKKY)  function  and  for  small  values  of  Rat)  the  (RKKY) 
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Table  I. 

Rare 


Earth 

S 

L 

a 

Tb 

3 

3 

-1/6 

Dy 

5/2 

5 

-1/15 

Ho 

2 

6 

-1/44 

Er 

3/2 

6 

1/33 

Tm 

1 

5 

1/6 

function  represents  a  very  poor  approximation  to  iso¬ 
tropic  interaction. 

In  order  to  obtain  an  estimate  of.  the  effect  of 
the  anisotropy  on  the  crystalline  exchange  energy  we 
investigated  the  energy  as  a  function  of  moment  orienta¬ 
tion  for  the  ferromagnetic  phase.  Using  an  IBM  7040  we 

ions 

summed  b  in  Eq.  (2)  over  all  of  a  hep  lattice  inside  the 
sphere  centered  at  a  and  of  radius  equal  to  50  nearest 
neighbor  spacings.  Prom  the  form  of  the  anisotropy  in 
Eq.  (2)  and  from  the  symmetry  of  a  hep  lattice  the  ferro¬ 
magnetic  energy  depends  only  on  the  angle  between  the 
moments  and  the  c-axis.  The  results  of  this  calculation 
is  given  in  Table  II  for  the  two  extreme  cases  in  which 
the  moments  are  along  the  c-axis  or  in  the  hexagonal 
plane . 

In  Table  II  it  can  be  seen  that  there  is  about  a 
10$  difference  in  energy  for  the  two  orientations,  and 
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Rare 

Earth 

Table 

Exchange 

II. 

Energy 

Along  c-axis 

in  Hex.  Plane 

Tb 

-O.I56 

-0.175 

Dy 

-0.106 

-0.121 

Ho 

-0.0704 

-O.O75O 

Er 

-0.0429 

-0.0395 

Tm 

-0.0213 

-O.OI55 

in  each  case  the  anisotropy  favors  the  orientation  which 
is  observed  experimentally  with  regard  to  whether  the 
moments  are  aligned  along  the  c-axis  or  in  the  hexagonal 
plane.  Although  there  are  other  effects  contributing  to 
this  such  as  crystal  field  terms,  on  the  basis  of  the 
exchange  energy  alone  one  can  explain  the  preferred  direc¬ 
tion  of  moment  orientation. 

Using  a  free  electron  gas  model  we  have  calculated 
exactly  the  Isotropic  and  first  order  non-isotropic  terms 
of  the  interaction,  and  on  the  basis  of  a  simple  ferro¬ 
magnetic  order  pattern  we  have  shown  that  the  anisotropy 
contributes  significantly.  This  calculation  by  no  means 
represents  the  true  physical  picture  since  plane  waves 
were  used  for  the  conduction  electrons.  Until  an  accurate 
treatment  is  done  with  the  appropriate  band  functions  the 
nature  of  the  interaction  will  remain  in  doubt. 
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ABSTRACT 

Ttfo  atomic  percent  of  the  various  tri valent  rare  earths  were  dis¬ 
solved  in  cerium,  and  the  electrical  resistivities  of  these  dilute 
binary  alleys  evaluated  in  the  range  1.5  to  29T°K,  in  order  to  determine 
whether  the  various  solute  magnetic  properties  might  be  reflected  in 
resistivity  measurements.  It  was  found  that  the  y  and  p  phases  at  room 
temperature  showed  a  roughly  linear  dependence  of  resistivities  on 
atomic  number  of  solute  species,  with  a  slight  effect  attributed  to 
solute-solvent  atomic  size  mismatches.  At  4°K  the  resistivities  of  pre¬ 
dominantly  a,  phases  were  found  to  be  influenced  by  the  presence  of  small 
amounts  of  p,  and  predominantly  p  phases  by  small  amounts  of  a.  Spin- 
disorder  resistivities  for  p  phase  were  found  to  be  a  nearly  linear 
function  of  atomic  nisnber  of  solute,  and  did  not  clearly  reflect  the 
influence  of  any  solute  magnetic  parameters.  A  slight  decrease  in  spin- 
disorder  resistivities  was  noted  in  the  direction  of  the  heavier  rare 
earth  solutes. 


INTRODUCTION 

The  element  cerium  can  exist  in  three  phases  between  room  tempera¬ 
ture  and  0°K.  After  being  annealed  at  high  temperature,  the  room  temper¬ 
ature  crystal  structure  is  fee  (y  phase),  with  a  lattice  constant  of 
5.16  A.  Cerium  in  this  form  has  about  0.9  4f  electron  per  atom.  On 


Work  performed  under  the  auspices  of  the  V.  S.  Atomic  Energy  Commission. 
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cooling,  the  7  phase  transforms  partially  to  th-  hep  form,  but  the 
number  of  4f  electrons  is  unchanged.  With  further  cooling,  the  untrans¬ 
formed  7  phase  undergoes  an  electronic  transformation  to  a  more  dense 
fee  form,  designated  c:  phase.  This  structure  has  a  lattice  parameter 

O 

of  4.85  A,  and  about  0.4  4f  electron  per  atom.  With  still  further  cool¬ 
ing,  the  0  phase  transforms  to  the  a  phase,  thus  undergoing  both  a 
crystallographic  and  an  electronic  transformation.  These  low-temperature 
cerium  phases  have  been  more  completely  described  by  Gschneidner.  ^ 

Repeated  thermal  cycling  of  annealed  samples  between  roam  and 

liquid  helium  temperatures  can  increase  the  amount  of  3  to  the  point 

■  (2) 

where  it  is  predominant  in  this  temperature  range.'  1  The  3  phase  of 
cerium  transforms  from  the  paramagnetic  to  the  antiferromagnetic  form, 
at  12.5°K.^ 

Previous  experimental  findings  on  cerium  containing  small  amounts 

of  other  rare  earths  suggest  that  the  electrical  conductivity  of  such 

binary  solid  solution  alloys  may  be  a  function  of  solute  species. 

(4) 

Gschneidner  et  al.  '  reported  that  the  temperature  of  the  7  -*  a  cerium 
transformation, which  involves  a  change  in  electronic  configuration  is 
strongly  dependent  on  the  specific  rare  earth  alloying  element  at  the 
2  a/o  level.  In  explaining  this  effect,  they  postulated  that  the  spin 
of  the  unpaired  4f  solute  atoms  was  primarily  responsible  for  the 
observed  dependence.  Their  conclusion  resulted  from  the  observation 
that  the  change  in  transformation  temperature  as  a  function  of  the  atonic 
number  of  the  rare  earth  solute  was  reciprocally  related  to  the  spin  of 
the  solute  atons.  The  mechanism  by  which  transition  temperatures  were 
lowered  was  believed  to  involve  an  indirect  interaction  between  the 
unpaired  4f  electron  of  the  cerium  and  those  of  the  solute  atons  via  the 

(4) 

conduction  electrons.'  '  This  suggestion  that  conduction  electrons  might 
play  a  part  in  the  observed  variation  of  transformation  temperature  as  a 
function  of  solute  species  in  cerium  led  to  the  present  experiment,  in 
which  the  resistivity  of  cerium  plus  2  a/o  rare  earth  additions  has  been 
determined  as  a  function  of  temperature  between  1.5  and  297°K.  The  pre¬ 
sent  paper  deals  with  the  resistivities  of  these  dilute  alloys  at  room 


■X- 

temperature  and  at  liquid  helium  temperatures.  Data  are  reported  as  a 
function  of  atonic  number  of  solute,  in  order  to  determine  the  possible 
effect  of  various  so.lute  magnetic  parameters  on  the  resistivities  of  the 
alloys.  Phases  considered  are:  y  at  297°K>  6  at  297°K>  3  at  4°K,  and 
a  at  4°K.  Only  the  3  alloys  at  297°K  were  thought  to  be  single  phase 
(i.e.,  phase  pure).  The  other  phases  listed  were  the  predominant  ones 
in  two-phase  mixtures. 

It  is  worth  emphasizing  that  the  cerium  alloys  dealt  with  here  were 
quite  dilute.  Only  one  solute  atom  was  present  for  every  50  solvent 
atoms. 

Only  one  other  investigation  of  the  effect  of  a  series  of  dilute 
solutions  of  rare  earths  in  another  rare  earth  on  the  electrical  resis¬ 
tivities  of  the  alloys  has  come  to  our  attention.  Mackintosh  and 

(5) 

Smidtv  /  have  reported  residual  resistivities  of  lutetiun  containing 

« 

small  percentages  of  gadolinium,  terbium,  dysprosium,  hoi  mi  urn,  and  erbium. 
•They  concluded  that  the  change  in  residual  resistivity  per  atomic  percent 
of  solute  added  is  a  function  of  the  factor  (g-l)  J(J+1),  where  g  is  the 
Lande  factor  and  J  is  the  localized  ionic  moment  of  the  solute. 

EXPERIMENTAL  PROCEDURE 

The  eleven  alleys  discussed  here  were  made  from  cerium  stock  which 
was  estimated  by  difference  to  be  99«75  w/o  pure.  The  major  metallic 
Impurities  were  0.1  w/o  lanthanum,  0.05  w/o  manganese,  and  0.05  w/o 
iron.  Major  non -metallic  impurities  were  255  PI®  carbon,  115  ppm  hydro¬ 
gen,  50  ppm  nitrogen,  and  80  ppm  oxygen.  The  resistivity  of  the  cerium 
stock  at  5«96°K  (the  normal  boiling  point  of  liquid  helium  at  Los  Alamos) 
was  7.14  pfl-cm. 

Most  of  the  alloys  were  arc-cast,  swaged  from  about  0.512  to  0.265 
in.  diameter,  heat  treated,  and  then  machined  into  the  final  rod-shaped 


* 

The  effects  of  the  normally  trivalent  rare  earth  solutes  are  reported 
here .  The  normally  divalent  elements  europium  ?nd  ytterbium  will  be 
discussed  at  a  later  date. 
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form.  Pinal  sample  dimensions  were  0.250  in.  diameter  by  1.125  in* 
long.  Potential  probe  separations  for  electrical  resistivity  measure¬ 
ments  were  about  1.5  cm. 

Arc-cast,  alloys  containing  dysprosium,  erbium  and  thulium  were  not 
swaged,  but  were  induction  cast  to  the  approximate  final  shape  in  a 
helium  atmosphere.  The  casting  technique  was  to  force  the  molten  alloy 
up  into  a  0.275  in.  ID  tantalum  tube  immersed  in  the  melt,  by  means  of 
helium  admitted  to  the  previously  evacuated  furnace  (man ©metric  casting 
method) .  The  tantalum  tube  was  subsequently  machined  from  the  casting, 
and  the  specimen  reduced  to  its  final  dimensions  by  further  machining 
after  heat  treatment , 

The  alloys  were  homogenized  and  annealed  in  the  following  manner. 
Each  specimen  was  wrapped  in  tantalum  foil  and  sealed  within  a  Pyrex 
capsule  under  0.5  atm  of  argon.  The  capsules  were  then  placed  in  a 
furnace  at  525  °C  (about  0.7  of  the  melting  point)  for  approximately  300 
hr.  At  the  end  of  this  period,  the  capsules  were  removed  from  the  fur¬ 
nace  and  cooled  rapidly  (about  100°K/min)  to  room  temperature.  The 
specimens  were  then  broken  out  of  the  capsules  and  machined  to  their 
final  dimensions.  Specimen  densities  were  determined  at  this  time,  and 
these  7-phase  values  are  reported  in  Table  I. 

The  machined  specimens,  four  at  a  time,  were  loaded  into  the  copper 
block  specimen  holder  of  the  resistivity  apparatus,  and  their  resistivi¬ 
ties  measured  at  297°K.  These  data  are  reported  as  the  gamma  curve  of 
Fig.  1(a). 

Cooling  to  liquid  helium  temperatures  was  then  carried  out  in  the 
following  manner.  The  assembly  containing  the  specimens  was  pre-cooled 
rapidly  in  liquid  nitrogen  and,  after  the  nitrogen  had  been  removed,  was 
quenched  with  liquid  helium.  The  liquid  nitrogen  pre-cool  was  utilized 
to  minimize  subsequent  liquid  helium  boil-off.  The  time  required  for 
cooling  from  297  to  4°K  was  typically  2  hr.  This  rapid  cooling  was 
utilized  to  favor  the  formation  of  a.  It  has  been  reported  that  slow 
cooling  enhances  the  formation  of  (3  phase,'  1  which  is  slower  than  7  to 
transform  to  a  on  further  cooling . ^ 
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Resistivity  data  taken  at  4°K  on  such  predominantly  alpha-phase 
samples  are  shown  by  the  "alpha"  curve  of  Fig.  1(b).  Pumping  was  then 
employed  to  reduce  the  temperature  of  the  liquid  helium  to  approximately 
1.5°K.  A  resistance  heater  imbedded  in  the  copper  specimen  holder  was 
used  to  warm  the  specimens  from  this  temperature  to  297°K  at  about 
O.^K/min. 

Specimens  to  be  evaluated  in  the  p -phase  condition  were  given  the 
following  t hemal  cycling  treatment  after  being  heat  treated  and  before 
being  inserted  in  the  resistivity  apparatus .  (  Each  was  sealed  within  a 
Pyrex  capsule  under  1  atm  of  helium  and  then  thermally  cycled  2k  times 
from  room  -  to  liquid  helium  -  to  roan  temperature .  The  procedure  for 
one  cycle  involved  immersing  the  capsule  in  liquid  helium,  holding  it 
in  the  liquid  for  15  min,  and  then  removing  it  for  warming  to  room  tem¬ 
perature  over  a  15  min  period.  Specimen  densities  were  determined  after 
2k  cycles,  and  these  p-phase  values  are  given  in  Table  I. 

The  p  samples  prepared  in  the  above  manner  were  leaded  into  the 
apparatus  and  their  resistivities  measured  at  297*K»  Brese  values  are 
reported  in  the  beta  curve  of  Fig.  1(a).  They  were  then  slowly  cooled 
at  a  rate  of  about  0. 5°K/min  to  about  40*K,  quenched  with  liquid  helium, 
and  the  helium  puntped  to  approximately  1.5*K.  Specimen  resistivities 
were  determined  at  the  liquid  helium  isotherm  of  h*K  and  arc  give®  by 
the  "beta"  curve  of  Fig.  1(b) .  As  was  the  caae  for  the  alpha  specimens, 
the  apparatus  beater  was  utilized  to  vans  the  specimens  to  297*K  at  a 
rate  of  about  ^.5cK/xain. 

Resistivity  data  were  taken  manually  during  all  runs,  a  potentio- 
(7) 

metric  circuit'1'  being  utilized.  Voltage  drops,  which  were  later  con¬ 
verted  to  resistivities,  were  read  with  a  Rubicon  Type  B  potentiometer . 

The  sensitivity  of  this  circuit  was  such  that  resistivity  changes  a n 
-8 

small  as  5*10 "  ohms  were  reliably  detected. 

Two  methods  of  measuring  temperature  were  used.  In  the  range 
between  room  temperature  and  l6°K,  temperatures  were  obtained  by  averag¬ 
ing  the  results  from  two  calibrated  copper-constaatan  theimocouples 
whose  measuring  junctions  were  b altered  to  the  copper  specimen  holder. 
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The  accuracy  was  estimated  to  be  ±  0.5°K.  In  the  range  between  50  and 
1.5°K,  a  calibrated  germanium  resistance  thermometer  was  used,  which  had 
an  estimated  accuracy  of  ±  0-l°K.  Its  temperature  measuring  element  was 
clamped  to  the  specimen  holder. 

RESULTS  AND  DISCUSSION 

Measured  resistivities  of  cerium-rich  alloys  at  297°K  as  a  function 
of  rare  earth  solute  species  are  shown  in  Fig.  1(a).  The  lower  curve 
represents  samples  in  the  annealed  condition.  Density  measurements, 
shown  in  Table  I,  indicated  that  most  of  these  specimens  were  primarily 
7  phase,  although  the  observed  trend  toward  increasing  density  with 
increasing  atonic  number  of  the  solute  was  not  continued  for  the  solutes 
erbium,  thulium  and  lutet-ium.  This  failure  to  follow  the  trend  was 
attributed  to  the  likely  formation  of  a  small  amount  of  beta  phase  in 
these  three  alloys.  The  upper  curve  of  Fig.  1(a)  represents  resistivi¬ 
ties  of  the  cerium  alloy  samples  at  297  °K  after  they  had  been  thermally 
cycled  24  times  to  4°IC  to  transform  them  as  completely  as  possible  to  |3. 
The  trend  in  density  measurements  (Table  I)  and  subsequent  resistivity- 
temperature  curves  indicated  that  the  cycled  specimens  consisted  entirely 
of  (3  at  room  temperature.  It  may  be  seen  in  Fig.  1(a)  that  the  resis¬ 
tivities  of  both  7  and  ^  phases  are  linear  functions  of  the  atomic  num¬ 
bers  of  rare  earth  solutes. 

The  results  shown  in  Fig.  l(a)  indicate  that  at  room  temperature  the 
magnetic  properties  of  the  various  rare-earth  solute  species  dissolved 
in  cerium  to  the  extent  cf  2  a/o  have  no  appreciable  effect  on  the 
electrical  resistivity  of  the  alloys.  The  slight  linear  rise  in  resis¬ 
tivities  as  a  function  of  atomic  number  of  solute  is  probably  attributable 
to  lattice  strains  resulting  from  differences  between  the  atomic  sizes 
of  the  solvent  and  solutes.  Since  in  both  the  7  and  (3  phases  the  cerium 
atoms  are  tri' valent, ^  the  effect  of  the  lanthanide  contraction  in  the 
other  trivalent  rare  earths  is  to  cause  the  greater  size  mismatch  to 
occur  with  the  heavier  rare  earth  solutes.  This  would  then  lead  to 
higher  resistivities  for  those  alloys  near  the  lutetium  end  of  the  series, 
as  was  found  experimentally. 
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Measured  resistivities  at  4°K  are  shown  in  Fig.  1(b).  The  lower 
curve  represents  measurements  on  specimens  cooled  rapidly  to  liquid 
helium  temperature  from  the  annealed  state.  It  might  be  inferred  from 
this  curve  that  the  resistivities  of  these  samples,  which  were  predomi¬ 
nantly  a  phase,  is  a  function  of  some  solute  property  which  relates  to 
a  minimum  near  neodymium  or  samarium.  However,  examination  of  the 
curves  of  resistivity  as  a  function  of  temperature  for  the  various  sam¬ 
ples  heated  slowly  from  4°K  indicated  that  this  was  not  the  case. 

Instead,  it  appeared  that  the  resistivity  variation  as  a  function  of 
atomic  number  of  solute  was  due  to  the  different  percentages  of  p  phase 
present  in  these  otherwise  a  alloys  at  4°K.  It  was  possible  from  exam¬ 
ination  of  the  heating  curves  to  make  a  crude  estimate  of  the  amount  of 
P  present  in  the  a  phase,  by  observing  the  magnitude  cf  the  resistivity 
change  at  the  Neel  point  of  the  £3  phase.  The  smallest  discontinuities 
occurred  at  about  samarium  or  gadolinium  in  the  series,  and  the  largest 
at  lanthanum  and  lutetium.  Thus,  with  regard  to  the  “alpha"  curve  in 
Fig.  1(b),  we  conclude  tht-t  A>ie  various  amounts  of  high-resistivity  p 
phase  present  in  the  predominantly  a -phase  specimens  may  have  obscured 
any  possible  contributions  due  to  magnetic  properties  of  the  rare  earth 
solutes . 

The  upper  curve  of  Fig.  1(b)  represents  resistivities  of  the  ceriusn 
alloy  samples  at  k°K  after  24  thermal  cycles  to  that  low  temperature. 

The  departure  of  this  curve  from  a  straight  line  can  largely  be  attri¬ 
buted  to  a  lack  of  phase  purity.  It  was  found  from  examination  of  the 
resistivity  vs.  temperature  curves  of  these  predominantly  p  samples  that 
some  a  phase  formed  on  cooling  at  lew  temperatures  but  reverted  to  p 
phase  on  warning,  so  that  the  samples  at  room  temperatures  were  always 
in  the  p  form.  The  amount  of  a  present  in  these  samples  at  4°K  (as 
•  deduced  from  the  resistivity-temperature  curves)  was  greatest  at  samarium 
and  least  at  lanthanum  and  lutetium.  Thus  the  presence  cf  some  low- 
resistivity  a  phase  in  the  predominantly  p  samples  would,  at  least  in 
part,  account  for  the  shape  of  the  curve. 

We  think  that  the  resistivities  of  the  predominantly  p  samples  at 
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4°K  can  be  roughly  corrected  for  c: -phase  effects,  since  the  percentages 
of  o  from  sample  to  sample  appeared  to  be  rather  small.  The  corrections 
would  tend  to  straighten  the  "beta"  line  shown,  so  that  a  straight  line 
would  approximate  the  corrected  data.  A  straight  line  suggests  that  at 
low  temperatures  the  magnetic  properties  of  the  rare-earth  solutes  in  |3 
cerium  have  no  major  effect  on  the  electrical  resistivities  of  the  alleys. 
A  rise  in  resistivities  toward  the  heavy  end  of  the  series  is  observed, 
as  was  seen  at  higher  temperatures,  and  may  be  attributable  to  size 
effects.  The  magnitude  of  the  effect  at  low  temperatures,  however,  is 
larger  than  was  found  at  room  temperature  (Fig.  1(a)). 

We  have  thus  far  discussed  only  the  measured  values  of  electrical 
resistivity  as  a  function  of  solute  species  in  our  attempt  to  determine 
the  effect  of  the  magnetic  properties  of  the  solute  atems  on  resistivi¬ 
ties  of  the  alloys.  Of  more  fundamental  significance  are  the  results  of 
an  examination  of  spin-disord2r  resistivity,  since  this  parameter  is  not 
as  subject  to  secondary  influences  as  is  measured  resistivity.  Spin- 
disorder  resistivity  is  defined  here  as  the  difference  between  the  para¬ 
magnetic  resistivity  and  the  residual  resistivity,  both  extrapolated 
graphically  to  0°K.  The  paramagnetic  resistivities  were  extrapolated 
from  the  linear  paramagnetic  regions  in  the  resistivity-temperature 
curves-  the  residual  resistivities  were  extrapolated  from  the  resistivity 
data  obtained  at  the  lowest  temperatures . 

Spin-disorder  resistivities  were  thus  obtained  from  the  resistivity- 
temperature  plots  of  the  predominantly  (3  samples,  which  ordered  magnet¬ 
ically  at  ~  17°K.  The  parameters  used  to  calculate  spin-disorder  resis¬ 
tivity,  and  the  results  obtained,  are  shown  in  Table  II.  Spin- 
disorder  resistivity  as  a  function  of’  atomic  number  of  solute  is  shown 
in  Fig.  2.  The  dip  in  the  curve  at  samarium  is  partially  due  to  the  pre¬ 
sence  of  a  relatively  high  percentage  of  a  phase,  as  discussed  earlier. 
Our  experimental  evidence  indicates  that  the  spin-disorder  resistivity, 
if  present  at  all,  is  lower  in  a  phase  than  in  {3.  However,  a  rough 
correction  to  remove  the  effects  of  the  various  amounts  of  alpha  present 
in  the  (3  specimens  did  not  change  the  overall  shape  of  the  curve.  It 
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thus  appears  that  the  true  curve  for  spin -disorder  resistivity  as  a 
function  of  atomic  number  of  solute  would  have  slight  dips  near  the  cen¬ 
ter  and  at  the  ends  of  the  series,  s omevhat  reminiscent  of  the  variation 
of  the  magnetic  parameter  L  (L+l)  across  the  rare  earth  series,  where  L 
is  the  total  orbital  angular  momentum.  The  magnitude  of  the  dips,  how¬ 
ever,  is  about  the  same  as  that  due  to  experimental  scatter.  The  only 
trend  clearly  established  in  this  experiment  was  a  slightly  negative 
slope  toward  the  heavier  rare  earths. 

The  present  findings  show  no  clear  variation  of  resistivities  with 
any  magnetic  parameter  across  the  rare  earth  series.  If  such  a  depend¬ 
ence  does  exist,  it  is  likely  that  it  could  be  unambiguously  detected  by 
the  present  technique  only  if  a  higher  percentage  of  solute  were  present. 

Such  solute  concentrations,  however,  were  not  required  in  the  experiments 

(5) 

of  Mackintosh  and  Smidt,'  who  added  small  percentages  of  heavy  rare 
earths  to  lutetium  and  found  significant  changes  in  residual  resistivi- 

p 

ties  which  they  correlated  with 'the  magnetic  parameter  (g-l)  J(J+1). 

(g) 

Smidt  and  Daane,  when  they  alloyed  rare  earths  with  other  rare  earths, 
found  large  variations  in  spin-disorder  resistivities  and  a  correlation 
with  the  term  (g-l)2J(J+l) .  Their  experiments  are  not  directly  compar¬ 
able  with  ours,  however,  since  they  examined  much  higher  solute  concen¬ 
trations.  They  also  found  that  the  magnetic  ordering  temperatures  varied 
greatly  with  alloying  content.  In  the  present  experiment  negligible 
variations  in  p -phase  Neel  temperatures  across  the  rare-earth  solute 
series  were  found,  an  observation  consistent  with  the  absence  of  any 
detectable  dependence  of  resistivity  on  magnetic  parameters. 

The  same  cerium-rich  compositions  discussed  here  were  found  by 
Gechneidner  et  alv  7  to  show  a  solute  spin-dependent  effect  on  the  y  -»  a 
transformation  temperatures .  It  is  not  clear  why  a  similar  effect  was 
not  observed  in  the  present  electrical  resistivity  data  if  the  conduction 
electrons  are  involved  in  the  coupling  of  the  4f  electrons. 

The  negative  slope  seen  in  the  plot  of  spin-disorder  resistivity  as 
a  function  of  increasing  solute  atomic  nucfoer  (Fig.  2)  is  rather  minor. 
It  may  be  related  to  the  observation^ 10 ^  that  the  energetically  deep  4f 
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shells  in  the  heavier  rare  earths  are  highly  localized  and  do  not  inter¬ 
act  strongly,  whereas  the  energy  levels  for  the  4f  electrons  in  the 
lighter  rare  earths  are  closer  to  those  for  the  5&  %  states,  and  thus 
may  interact  to  some  extent.  From  this  view,  the  magnetic  interaction 
of  solvent  and  solute  atoms  by  the  Friedel-Rocher  mechanism  of  virtual 
bound  states (-^,14)  be  expected  to  be  less  with  heavy  rare  earth 

solutes,  and  thus  decrease  spin-disorder  resistivity  when  these  elements 
are  present  in  cerium.  It  is  also  possible  that  the  negative  slope  in 
Fig.  2  may  be  due  to  changes  in  band  structure  or  Fermi  level  as  alloy¬ 
ing  elements  are  added  to  cerium. 

SUMMARY 

Two  atomic  percent  of  the  various  tri valent  rare  earths  were  dis¬ 
solved  in  cerium,  and  the  electrical  resistivities  of  these  dilute 
binary  allqyw  detemined  in  the  range  1.5  to  297°K.  Resistivities  are 
reported  for:  predominantly  y  phase  at  297°K,  3  phase  at  297 °K,  pre¬ 
dominantly  a  phase  at  4°K,  and  predominantly  (3  phase  at  4°K. 

Resistivities  at  297°K  for  both  y  and  3  phases  were  found  to  be  a 
linear  function  of  atomic  number  of  solute,  with  a  slight  increase 
toward  the  heavy  rare  earths  being  attributed  to  atomic  size  effects. 

Resistivities  for  predominantly  a-phase  samples  at  4°K  were  found 
to  be  influenced  significantly  by  the  percentages  of  3  phase  retained  on 
cooling.  3- phase  resistivities  at  4°K  were  also  affected  by  small 
amounts  of  a  phase  present. 

Spin-disorder  resistivities  were  detemined  for  the  predominantly  3 
samples  and  were  found  not  to  reflect  the  influence  of  solute  magnetic 
parameters  that  had  been  suggested  by  the  work  of  others.  A  slight  down¬ 
ward  trend  in  spin-disorder  resistivities  was  detected  with  increasing 
atomic  number  of  solute,  and  is  considered  in  terns  of  a  virtual  bound 
state  model. 
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Table  I.  Room- temperature  gemma  and  beta  phase  densities  of 

cerium  containing  2  a/o  of  various  other  rare  earths. 
(Density  determination  by  an  immersion  method.) 


Solute 

La 

Pr 

Nd 

Sm 

* 

7 -phase  density 

6.73 

6.75 

6.76 

6.76 

(g/cm3) 

p~phase  density 

6.67 

6.69 

6,70 

6.71 

(g/cm3) 

Density  diff. 

0.06 

0.06 

0.06 

0.05 

Gd  lb  Dy  Ho  Er  lin  Lu 

6.77  6.77  6.78  6.78[6.77|[6.78][6.77] 
6.71  6.72  6.72  6.72  6.73  6.73  6.74 
0.06  0.05  0.06  0.06  o.o4  0.05  0.03 


Brackets  indicate  uncertainty  in  phase  purity. 


Table  II.  Resistivities  of  cerium,  predominantly  beta  phase,  containing 
2  a/o  of  various  other  rare  earths.  (Units  of  p  are  tiQ-cm.) 


Solute 

La 

Pr 

Nd 

Sm 

Gd 

n> 

Dy 

Ho 

Er 

Tn 

Lu 

p  at  4°K 

16.8 

15.0 

16.6 

19.0 

24.6 

25.0 

24.4 

27.5 

24.0 

26.2 

30.1 

Paramagnetic  p 

4-9.8 

47.7 

47.8 

46.9 

52.5 

54.0 

52.7 

56.3 

53.2 

55.5 

57-0 

extrap,  to  0°K 

Residual  p 

15.1 

12.9 

14.4 

17.0 

22.5 

22.9 

22.6 

25.6 

22.1 

24.4 

28.2 

extrap. to  0°K 

Spin-disorder  p 

34.7 

VJJ 

4=- 

OO 

33.4 

29.9 

29.8 

31.1 

30.1 

30.7 

31.1 

31.1 

28.8 

at  0°K 


* 


The  spin-disorder  resistivity  was  obtained  by  subtracting  the  residual 
resistivity  from  the  paramagnetic  resistivity,  both  extrapolated  to  0°K. 
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La  Ce  Pr  Nd  Pm  Sm  Eu  Gd  Tb  Dy  Ho  Er  Tm  Yb  Lu 
ATOMIC  NUMBER  OF  SOLUTE 


Fig.  Spin-disorder  resistivity  of  beta  cerium  (containing  2  a/o 
additions  of  various  other  rare  earths)  as  a  function  of  the 
atomic  nusfcer  of  solute.  Quotation  marks  around  the  word  beta 
indicate  m  uncertainty  in  the  phase  purity  of  these  alloys. 


- 
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SPIN  DEPENDENCE  OF  THE  ELECTRICAL 
RESISTIVITIES  OF  GADOLINIUM  ALLOYS* 
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ABSTRACT 

The  electrical  resistivities  of  various  intra- rare -earth 
alloys  of  gadolinium  at  4.2°,  77°,  and  298°K  are  discussed  with  respect 
to  the  theory.  Two  spin  effects  are  shown  in  the  scattering  of  the 
conduction  electrons . 
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We  have  determined  the  resistivities  of  various  gadolinium 
alloys  at  4.2°,  77°>  and  298°K.  The  electrical  data  permit  a  critical 
study  of  the  scattering  of  the  conduction  electrons.  The  following 
conclusions  are  particularly  important  and  thought  to  be  valid  for  all 
intra-rare-earth  alloys  having  localized  magnetic  moments. 

(1)  The  increases  in  the  residual  resistivity  of  Gd  are  strongly 
spin  dependent.  The  impurity  scattering  of  electrons  is  due  predomi¬ 
nantly  to  a  spin  inhomogeneity  effect. 

(2)  At  temperatures  above  the  Curie  point,  the  spin-disorder 
resistance  is  so  preponderant  that  the  impurity  effect  on  scattering 
becomes  obscure  and  the  intrinsic  effect  of  spin  disorder  dominates 
the  scattering. 

(3)  If  the  scattering  caused  by  the  two  spin  effects  just  stated 
can  indeed  be  suitably  treated  by  the  Born  approximation  method,  as  has 
been  done  by  DekJcer,  ^  the  spin-disorder  resistivity  of  Gd,  pg^,  above 
the  Curie  point  would  he  much  larger  than  the  value  obtained  by 


extrapolation . v 
present  theory. 


(2,3) 


Such  a  discrepancy  reflects  an  inadequacy  of  the 


We  chose  Gd  as  the  base  metal  for  its  simple  S^/0  ground 
state.  The  single  magnetic  transition  and  its  Curie  temperature 
(0  =*  289°K)  were  also  considered  to  be  desirable.  The  solute  elements 

V 

examined  include  La  and  all  rare  earths  except  Pm,  Eu,  and  Yb. 

At  4.2°  and  77°K,  all  solutes  were  observed  to  increase  the 
resistivity  of  Gd  according  to  the  Jfordheim  approximation.  By  plotting 
Ap/(l-c)  vs.  c,  vo  can  evaluate  the  atomic  resistivity  increase  caused 


~ -J  , 
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by  impurity  in  terms  of  4>/c(l-c),  where  «  Pgjj^y  “  Pq^  and  c  is 

the  atomic  concentration  of  the  solute.  In  Fig.  1,  the  values  of 

o  p 

4>/c(l-c)  deduced  from  the  data  at  4.2°K  are  plotted  against  S' ,  where 
S  has  the  value  of  the  difference  between  the  atomic  numbers  of  Gd  and 
the  solute  element.  It  is  seen  that  all  experimental  points  except 
that  of  Ce  fall  upon  two  straight  lines  parallel  to  each  other  within 
the  experimental  error.  Hence  ve  obtain  an  empirical  formula  of  the 
form 

4>  »  c(l-c)(A  +  BS2)  (1) 

with  apparently  different  A. values  but  the  same  B  value  for  the  light 

and  heavy  lanthanides  c  Equation  (l)  is  also  valid  at  77°K,  "but  at 

298°K,  nearly  all  4)  disappeared  as  the  alloys  showed  resistivities 

within  +  2f>  of  the  value  (155  jifl-cm)  for  pure  Gd. 

The  results  at  298°K  point  out  the  difference  between  the  Gd 

(b) 

and  certain  noble  metal  alloys  7  and  imply  a  strong  spin  dependence 
of  the  scattering  mechanism  leading  to  4>*  The  former  conclusion  is 
also  supported  by  the  fact  that  all  rare  earth  ions  under  considera¬ 
tion  (except  Ce,  see  Fig.  l)  have  equivalent  electron  configurations 

t  Sj\ 

leading  to  a  pseudotrivalence .  '  Therefore  the  Gd  alloys  are 

expected  to  show  little  resistance  due  to  electrostatic  scattering,  as 

is  evidenced  by  the  small  Intercepts  of  the  oest-fit  lines  in  Fig.  1 

2 

and  BS  is  to  be  identified  with  the  resistivity  due  to  exchange 
scattering  with  S  becoming  the  difference  between  the  spins  of  Gd  and 
the  solute  ions. 


ivryr* 
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This  interpretation  of  Eq.  (l)  is  in  agreement  with  Dekker’s 
analysis  for  the  case  that  the  coupling  constant  G  of  the  indirect 
spin  Interaction  between  the  conduction  electrons  and  the  magnetic 
shells  is  invariant  with  the  spin  operator  of  the  magnetic  ions  S. 

The  Hamiltonian  of  the  interaction  is  written  as  -2G5(ri 
where  s^  is  the  spin  operator  of  a  conduction  electron  at  r^.  The 
question  as  to  whether  G  is  constant  or  not  has  been  discussed  by 
Liu^  ^  and  Watson  and  Freeman .  ^ 

Although  the  Indirect  spin  interaction  is  essential  in  the 
consideration  of  the  exchange  scattering,  the  interaction  per  se  does 
not  cause  the  scattering.  Thus  Kasuya^  has  predicted  a  zero  resis¬ 
tivity  for  perfectly  pure  Gd  at  0°K.  The  real  cause  of  the  exchange 
scattering  is  the  inhomogeneity  of  the  spin  lattice  produced  by  the 
solute  atoms  and  the  scattering  occurs  only  at  the  solute  sites  where 

the  perturbations  are  centered.  These  points  have  been  overlooked  by 
(3) 

some  authors.' 

p 

It  should  also  be  cited  that  the  simple  relation  ~  S 

has  never  been  observed  in  the  alloys  of  Fe,  Co,  and  Ni.  This  is 

probably  because  the  wave  functions  of  the  3d  electrons  are  not  so 

localized  as  those  of  the  kt  electrons.  Also,  in  the  3d  transition 

metal  alloys,  the  Coulomb  scattering  potentials  are  not  necessarily 

small  and  complications  arise  such  as  the  resonance  scattering  caused 

(9) 

by  the  virtual  states.  ' 

To  illustrate  the  difficulty  encountered  in  the  reconcilia¬ 
tion  of  the  observed  in  the  Gd  alloys  with  the  spin-disorder 
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resistivity,  ve  recall  that  both  Dekker  and  de  Gennes  and  Friedel^0^ 
have  arrived  at  a  S'(S*  +  1)  dependence  of  pgd  for  Gd,  where  S‘  is 
the  atomic  spin  of  Gd.  Further,  from  Dekker* s  analysis,  an 
expression  may  be  obtained  comparing  p  with  the  resistivity  increase 
caused  by  the  spin  inhomogeneity  effect,  p^,  in  a  Gd  alloy  at  0°K 


P 


sa<T  »  ec> 


S*(3»  +  1) 


cS 


(2) 


Since  Ps^(0)  Bay  "be  estimated  from  the  results  shown  in  Fig.  1,  pgd 

can  be  calculated  from  Eq.  (2).  (Diking  Gd^^La^  as  a  typical  example, 

Zip  at  4.2°K  is  11  pH- cm  and  pgi(0)  is  approximately  10  pa-cm.  With 

S*  ■  7/2,  pBd  is  calculated  to  be  430  pa-cm,  as  compared  with  the 

2  3 

extrapolated  value  of  111  -120  pa-cm. 

This  large  discrepancy  between  the  calculated  and  extrapo¬ 
lated  values  of  pgd  demands  a  careful  re-examination  of  the  theory. 

One  possible  explanation  might  be  that  the  scattering  involving  spin 
flip  is  merely  a  mathematical  consequence  induced  by  the  spin  step-up 
and  step-down  operators  components,  ~-(s  S  +  s  S.),  of  the  s<S 

d.  +  -  -  +  —  — 

coupling.  In  reality,  the  paramagnetic  state  may  show  short-range 
order  of  the  spins,  thereby  allowing  only  the  elastic  scattering  due  to 
the  Bz^z  component  of  the  coupling  in  the  temperature  range  from  which 
the  extrapolated  value  of  pgd  has  been  obtained.  If  this  is  true,  pgd 

p 

would  depend  upon  S*  ,  instead  of  S*(3*  +  1).  By  assuming  equal 

p  p 

probability  for  all  possible  values  of  S'  ,  we  have  S*  «  21/4  and 
the  calculated  value  of  pgd  is  143  na-cm.  !Ehe  discrepancy  would  then 


» 


r.  t  "C 


“«5S3E222 


3T 


48 


be  reduced  to  a  satisfactory  level.  To  test  this  possibility, 
resistance  measurements  vere  conducted  on  a  Gd  sample  at  elevated 
temperatures.  The  resistivity  of  Gd  was  observed  to  increase  para- 
bolically  vith  temperature  from  500°K  to  a  value  of  190  microhm-cs  at 
1050°K,  and  remain  almost  unchanged  between  1050°  aDd  1250°K.  There 
was  no  indication,  at  least  up  to  1250°K,  that  the  resistivity  of  Gd 
would  rise  sharply  to  reach  the  predicted  value  of  pg^  Including  spin 
flip. 


The  author  is  indebted  to  Dr.  P.G.  KLemens  for  his  valuable 
contribution  to  the  discovery  of  the  discrepancy  in  the  calculation 
of  pgd  and  comments  on  the  manuscript  and  to  Dr.  W.J.  Carr,  Jr,  for 
stimulating  discussion. 
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Fig.  1— A  linear  relationship  is  seen  between  4p/c(l-c)  and  S 
at  4. 2°K  for  both  light  (La-Sm)  and  heavy  (Tb-Lu)  rare  earth 
solutes  in  Gd.  The  dashed  lines  indicate  a  possible  explanation 
of  the  large  deviation  from  the  linear  relationship  of  the  Ce 
point.  Instead  of  6,  S  is  approximately  6. 6,  corresponding 
to  60%  of  the  Ce  ions  in  the  tetravalerit  state. 
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MAGNETIC  TRANSFORMATION  IN  HEAVY  RARE-EARTH  ALLOYS 

WITH  EACH  OTHER 


R.  M.  Bozorth  and  R.  J.  Gambino 


IBM  Watson  Research  Center 
Yorktown  Heights,  New  York 


ABSTRACT 

We  have  measured  the  Curie  (T^,)  and  Neel  (T^)  points  of  the 
binary  intra-rare-earth  alloys  Gd-Lu,  Dy-Ho,  Dy-Tm,  Er-Tm,  and 
Tm-Lu,  in  addition  to  these  previously  reported  (Atlantic  City, 
Nottingham).  There  are  maxima  in  the  Curie  point  vs  composition 
curves  for  Ho-Er,  Ho-Tm,  and  Er-Tm,  minima  for  Dy-Tm  and 
Ho-Tm,  when  the  average  concentration  of  4f  electrons  is  near  10.  5. 
The  Tm-Lu  alloys  are  ferromagnetic  when  Tm-rich,  and  distilled 
T.n  has  a  Curie  point,  N&el  point,  and  intermediate  transformation 
point.  The  Lu-rich  alloys  are  all  non-ferromagnetic  in  low  fields 
and  some  become  ferromagnetic  at  a  critical  field,  when  they  have 
sharply  increased  susceptibility  with  onset  of  hysteresis.  Repre¬ 
sentative  data  for  various  alloy  systems  are  shown  to  deviate  only 
slightly  from  a  straight  line  when  plotted  against  the  E/3  power  of  the 
average  de  Cannes  factor,  as  discussed  by  Koehler. 
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METHODS 

Materials  were  prepared  by  arc  melting  and  annealing  in  closed 

1  2 

containers  and  then  tested  for  homogeneity  as  previously  described.  ’ 

o 

The  anneal  was  normally  1050  C  but  alloys  containing  Tm  were 

annealed  at  850°C„  The  distilled  Tm  was  prepared  by  sublimation  in 

a  vacuum  of  10  mm  at  1  500°C. 

Magnetic  measurements  were  made  with  the  pendulum  magne- 

3 

tometer  previously  described.  Although  several  methods  were  used 
for  determining  the  Curie  point  (T  ),  the  method  generally  employed 
depends  on  finding  the  point  of  inflection  of  the  moment  (O')  vs  tem¬ 
perature  (T)  curves  for  two  or  throe  relatively  low  fields  and  extrap¬ 
olating  these  points  to  zero  field.  Normally  the  point  of  inflection  is 
determined  visually  but  it  is  possibly  also  to  plot  ACT/ AT,  using  the 
AC  and  AT  between  neighboring  points  on  the  C,  T  curve,  and  so 
finding  the  temperature  of  maximum  slope.  This  "inflection"  method 
as  applied  to  Ho  and  Er  gives  the  same  values  of  T  as  those 
previously  reported  in  the  literature,  and  Dr.  R.  Joenk  has  shown 
that  the  method  is  consistent  with  molecular  field  theory.  In  the  case 

of  Tm  the  inflection  point  comes  at  25°K,  3°  above  the  temperature 

4 

previously  reported,  and  will  be  discussed  below. 

The  Neel  points  are  readily  detected  as  sharp  maxima  in  the 
C,  T  curves,  extrapolated  to  zero  field.  The  difference  in  the  char¬ 
acter  of  the  singularity  in  the  C,  T  curve  at  T^,  and  T^  is  especially 
apparent  in  the  Gd-rich  alloys,  in.  some  of  which  T^,  and  T^  coincide 
and  are  properly  considered  to  be  Curie  points,  whereas  other  alloys 
show  on  cooling  first  the  sharp  peak  of  a  Neel  point  and  at  lower 
temperatures  these  may  or  may  not  show  a  true  Curie  point. 

RESULTS 

In  Figure  1  we  plot  T^.  and  T^  against  the  average  number,  n, 
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of  4f  electrons  per  atom  of  alloy.  As  previously  noted  is  nearly 
a  universal  function  of  n  as  long  as  the  n's  of  the  component  elements 
are  not  too  different.  For  Gd-Lu,  however,  there  is  a  large  departure 

from  most  of  the  other  systems  such  as  Ho-Er.  Following  Koehler 

6  “  Z  /  3 

et  al.  this  departure  disappears  when  is  plotted  against  G  , 

where  G  is  the  average  of  the  de  Gennes^  factors: 


G  =  (g-i)  J(J+i) 


G  =  C1G1  +  C2G2  * 


Cj  and  c^  being  the  atomic  fractions  for  which  G  =  G^  and  G^  .  Fig¬ 
ure  2  show 2  such  a  plot  for  many  of  our  alloys  and  also  for  a  number 

7 

reported  by  the  group  at  Oak  Ridge  National  Laboratory.  The  de¬ 
parture  of  some  of  the  points  from  the  straight  line  appears  to  be  well 
outside  the  experimental  error;  recently  we  have  confirmed  the  deter¬ 
mination  of  T^  for  one  of  the  Gd-Er  alloys*  from  the  presumed  best 
line.  There  thus  appears  to  be  some  deviation  from  the  empirical 
2/3  power  relation  not  accounted  for  by  experimental  error  or  im¬ 
purities. 

Data  used  in  the  determination  of  the  Nfcel  and  Curie  points  of 
tw'o  Tm-Lu  alloys  are  shown  as  examples  in  Figure  3.  Tr,  for 

V 

O 

Tm  ,  Lu  .  ,  if  it  exists  at  all,  is  below  2  K,  but  evidence  of  ferro- 

U%  O  U*  t 

magnetism,  or  of  some  change  in  structure,  is  shown  by  the  change  in 
slope  of  the  magnetization  curves  of  Figure  4, 

The  values  of  T  plotted  in  Figure  i  are  distinctly  characteristic 

c 

of  each  alloy  system.  Greater  departures  from  linearity,  in  the  curves 
joining  two  elements,  are  found  in  the  middle  of  the  diagram  near,  the 

g 

concentration  of  4f  electrons  of  10.  5,  known  to  be  the  point  at  which 

the  4f  electron  cloud  changes  from  prolate  to  oblate  and  as  a  result 

9  10 

the  crystal  anisotropy  and  magnetostriction  change  in  sign.  The 
fractional  increase  in  T  of  the  alloy  as  compared  to  that  of  the  com- 

v 

O  o 

ponent  elements  is  high  in  the  Ho-Er  alloys(35  ,  cf.  20  )  and  even 
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higher  in  the  Er-Tm  alloys  (43°,  cf.  20°  and  25°),  On  the  contrary 
T_  is  depressed  by  alloying  Dy  and  Tm  and  T  is  then  lowered 
from  85°  and  25°  to  7°K. 

There  is  special  interest  in  the  element  Tm.  Previous  measure¬ 
ments,  based  on  the  existence  of  hysteresis,  indicated  that  T  =  22°. 

L/ 

Recently  we  have  measured  the  loops  of  material  obtained  from  the 
manufacturer  and  annealed,  and  the  same  material  as  distilled  and 
analyzed,  at  i.  25°K  with  =  25  kOe.  As  may  be  seen  in  Figure  5 

the  former  material  showed  hysteresis  with  a  coercive  force  of  about 
4  kOe  and  a  sharp  increase  in  0"  as  H  exceeded  20  kOe.  Spectro¬ 
scopic  analysis  showed  0.  6  at.  %  Ta  and  0.  4  at.  %  O.  On  the  contrary 
the  distilled  material,  c<  utaining  0.  02  at.  %  Ta  and  0.  4  at.  %  O, 
showed  no  such  increase  up  to  25  kOe.  Measurements  by  T.  R.  Mc¬ 
Guire  in  fields  to  50  kOe,  however,  showed  a  sharp  increase  in  O' 
beginning  at  about  30  kOe  and  still  rising  rapidly  at  50  kOe,  when  the 
moment  was  about  2.  5  Bohr  units  per  atom.  Apparently  the  4,  3,  4,  3 
structure  found  by  neutron  diffraction^  in  zero  field  is  broken  down 
in  high  fields.  (The  theoretical  gj  is  7.  0.  ) 

The  point  of  inflection  of  the  CT,  T  curve  for  Tm  is  found  to  be 
25°K,  as  extrapolated  to  H  =  0.  The  slope  of  the  curve  (Figure  6)  was 
determined  for  distilled  material  with  some  accuracy  by  measuring 
both  ACT  and  AT  between  neighboring  points  taken  for  H  =  2  kOe. 

A  plot  of  the  (negative)  slope  (Figure  7)  shows  a  well-defined  maxi¬ 
mum  at  26°K.  At  33°K  a  sharp  decrease  in  the  C,  T  curve  is  ob¬ 
served  such  that  the  slope  then  rises  suddenly  to  a  new  higher  maxi¬ 
mum,  indicative  of  a  departure  from  the  4,  3,  4,  3  arrangement  and 
the  formation  of  a  new  magnetic  structure.  Sharp  decreases  of  this 
kind  were  found  for  three  different  specimens  of  Im  of  different  origin. 

The  Curie  points  of  Tm-rich  Tm-Lu  alloys  decrease  rapidly  to 
zero,  and  are  not  ferromagnetic  above  a  concentration  of  about  30%  Lu. 


—  2/3 

Neel  points  were  determined  to  about  90%  Lu  and  follow  the  G  law 

with  some  accuracy  (see  Fjgure  2). 

R  eferences 

1.  R.  M.  Bozorth  and  R.  Gambino,  Proc.  Intern,  Gonf.  on  Mag¬ 
netism,  Nottingham  (1964),  in  press,  (p.  263). 

2.  R.  M.  Bor  *('i  and  J.  C.  Suits,  J.  Appl.  Phys.  35,  1039  (1964). 

3.  Bozorth,  Williams,  and  Walsh,  Phys.  Rev.  103,  572  (1956). 

4.  D.  D.  Davis  and  R.  M.  Bozorth,  Phys.  Rev.  118,  1543  (I960). 

5.  P.  G.  de  Gennes,  J.  phys.  radium  23,  510  (1962). 

6.  Koehler,  Wollan,  Child,  and  Cable,  Rare-Earth  Research, 

(K.  S.  Vorres,  Ed.  )  Gordon  and  Breach,  New  York,  p.  199 
(1964).  See  also  Weinstein,  Craig,  and  Wallace,  J.  Appl.  Phys. 
34,  1354  (1963),  for  the  use  of  the  3/2  power  law  of  vs  com¬ 
position  of  solid  solutions  with  yttrium. 

7.  Summarized  by  H.  R.  Child,  Thesis,  University  of  Tennessee 
(1965). 

8.  K.  W.  H.  Stevens,  Proc,  Phys.  Soc.  65,  209  (1952). 

9.  R.  >T.  Elliott,  Phys.  Rev.  1 24,  346  (.1961), 

10.  Tsuya,  Clark,  and  Bozorth,  Proc.  Intern.  Conf.  on  Magnetism, 
Nottingham  ( 1 964),  in  press. 

11.  Koehler,  Cable,  Wollan,  and  Wilkinson,  Phys.  Rev.  1 26,  1672 

(1962). 


56 


>; 

Q 

i 

XI 

t- 


o 

X 
l 

XI 

h- 


k. 

UJ 

I 

xs 

h- 


■f 

XI 

h* 


UJ 

I 

o 

X 


o 

X 

I 

>» 

Q 


I 

XX 

I- 


>"  w 
I  Q 
XX  _J 

X 

6  H 


to 

o 

rH 
r— H 

«J 

M 

(X) 

CJ 

•rH 

X 

a> 

£ 


bo 

•  rH 


Me‘3aniva3dW3i  133N  ao  3iano 


NUMBER  OF  4f  ELECTRONS 


for  Tm.  T„  for 


0  5  10  15  20  25 


H (KOE) 


Hysteresis  loops  for  Tm  specimens  as 
received  and  annealed,  and  as  distilled. 
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at  point  of  inflection,  and  discontinuity  at  33  K. 


63 


MAGNETIC  PROPERTIES  OF  RARE  EARTH- THORIUM  ALLOYS  ^ 

t  Research  sponsored  by  the  U.  S.  A.  E.  C.  under 
contract  with  the  Union  Carbide  Corporation. 

W.  C.  .Koehler,  H.  R.  Child  and  J.  W.  Cable 

Oak  Ridge  National  Laboratory 
Oak  Ridge,  Tennessee 

ABSTRACT 

Neutron  diffraction  studies  have  been  carried  out  on  a  number  of 
alloys  of  thorium  with  the  heavy  rare  earth  metals  Tb,  Ho,  and  Er  as 
part  of  an  investigation  of  the  nature  of  the  exchange  interaction  in  these 
materials.  At  room  temperature,  and  below,  the  heavy  rare  earths  in 
question  crystallize  in  the  simple  hexagonal  close  packed  structure; 
thorium,  in  the  cubic  close  packed  structure.  The  interatomic  distances 
in  thorium  are  intermediate  between  those  in  Dy  and  Ho,  and  relatively 
wide  regions  of  solid  solution  have  been  observed.  Alloys  containing 
more  than  about  80  a/o  rare  earth  are  found  in  the  hep  structure,  and 
those  with  less  than  about  50  a/o  are  found  in  the  ccp  arrangement. 

The  three  outer  electrons  of  the  rare  earth  metals  (5d  6s  )  are  pre¬ 
sumed  to  form  conduction  bands  and  it  is  generally  accepted  that  mag¬ 
netic  coupling  in  the  rare  earth  metals  occurs  via  an  indirect  interaction 
involving  a  polarization  of  the  conduction  electrons.  Thorium,  a  tetra- 
valent  diluent  (outer  electron  configuration  6d'7s  )  provides  a  means  of 
changing  the  effective  number  of  conduction  electrons  on  alloying  of  the 
rare  earths.  In  contrast  to  the  effect  of  trivalent  yttrium  and  lutetium 
as  diluents,  small  additions  of  thorium  tend  to  stabilize  the  ferromagnetic 
state  in  the  hep  structure.  For  example,  the  alloy  10Th-90Ho  transforms 
spontaneously  from  the  helical  to  the  ferromagnetic  state  and  in  the  alloy 
20Th-80Ho  the  helical  configuration  was  not  detected.  In  the  Th-rich 
region  alloys  20Tb-80Th,  and  30Tb-70Th  have  so  far  been  studied.  Both 
of  these  exhibit  diffraction  patterns  characteristic  of  short  range  anti¬ 
ferromagnetic  ordering  at  4.  2°K. 

Attempts  at  correlation  of  these  results  with  the  Ruderman-Kittel- 
Kasuya- Yosida  theory  have  been  made  but  with  indifferent  success. 
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PARAMAGNETIC  STUDIES  OF  HOLMIUM  BY  NEUTRON 
TOTAL  CROSS  SECTION  MEASUREMENTS 

Marieta  Mattos 

Institute)  de  Energia  AtSmica 
Sao  Paulo,  Brasil 

ABSTRACT 

The  total  neutron  cross  section  for  holmium  was  measured  for 
neutron  wavelengths  between  0.  20  and  9.  00  Angstroms.  The  IEA 
crystal  spectrometer  was  used  together  with  a  mechanical  velocity 
selector  to- eliminate  higher  order  reflections  from  the  crystal.  Powder 
samples  of  Hq203  were  prepared  in  high  purity,  taking  particular  care 
to  eliminate  contamination  by  the  rare  earths  of  high  neutron  cross 
section.  '  Although  paramagnetic  form  factors  are  usually  determined 
by  differential  cross  section  measurements  with  one  neutron  wave¬ 
length,  accurate  total  cross  section  measurements  for  a  wide  range  of 
known  wavelengths  give  definite  information  on  the  radius  of  the  4  i 
electron  orbits.  Our, measurements  indicate  a  radius  different  from 
the  one  used  by  Blume,  Freeman  and  Watson,  who  calculated  the  form 
factor  assuming  hydrogenic  wave  functions.  The  possible  choice  of 
different  screening  constant  will  be  discussed. 
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PARAMAGNETIC  STUDIES  OF  HOLMIUM  BY  NEUTRON  TOTAL 
CROSS  SECTION  MEASUREMENTS 

Marie ta  C*  Mattos 

Nuclear  Physios  Division,  Instituto  de  Energia  At&nica, 
Cidade  UniversitAria,  SSo  Paulo,  Brasil 


RESUMO 

Foi  medida  a  seog&o  de  choque  total  do  hftLmio,  para 
neutrons  de  coraprlmsntos  de  onda  entre  0,20  e  9,00  Angstroms  . 
Foi  usado  o  espectr6hietro  de  cristal  do  IEA  acoplado  a  urn  sele- 
tor  meoAhico  de  ve].ocidade  para  eliminar  reflexSes  de  ordem  su¬ 
perior  do  cristal.  Amostras  de  Ho^  em  p6  foram  preparadas  em 
alto  grau  de  pureza,  tomando-se  cuidado  especial  em  eliminar  con 
taminag&o  por  terras  raras  de  alta  secgSo  de  choque,  para  neu¬ 
trons.  Ainda  que  os  fatores  de  forma  psraraagndtica  sejam  usual- 
mente  determinados  por  medidas  de  secgffo  de  choque  diferencial 
para  neutrons  de  um  unico  comprimento  de  onda,  medidas  precisas 
da  secgfio  de  choque  total  num  grande  intervalo  de  comprimentos 
de  onda  conhecidos  dSo  informagffes  definidas  s$bre  o  raio  da  <5r- 
bita  eletrfinica  4f.  Nossas  medidas  indicam  um  raio  diferente  da 
quele  usado  por  Blume,  Freeman  e  Watson,  que  calcularam  o  fa- 
tor  de  forma  assumindo  fungCes  de  onda  hidrog&iioas.  SerA  dis- 
cutida  a  posslvel  escolha  de  diferentes  oonstantes  de  screening. 


La  section  efflcaoe  totale  de  l’holmium  a  AtA  mesurAe  pour 
les  neutrons  de  longueur  d’onde  entre  0,20  et  9/00  Angstroms.  Le 
spcctronAtre  A  cristal  de  1’lEA  a  At  A  employA  acouplA  A  un 
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s^lecteur  m^canique  de  vitesse  pour  dliminer  les  reflexions 
d’ordre  sup^rieur  du  cr  is  talc  Des  echantillons  de  H02O5  esipoudre 
ont  etd  prepares  avec  haute  puret£,  en  faisant  sp^cialement  atten 
tion  pour  dliminar  la  oontamination  pax  des  terres  rares  de 
haute  section  efficace  pour  les  neutrcns.  Bien  que  les  facteurs 
de  forme  soient  ordinalrement  determines  par  des  mesuras  de  sec¬ 
tion  effio&ce  differentielle  k  une  seule  longueur  d’onde  du  neu¬ 
tron,  des  mesures  precises  de  la  section  efficace  totale  dans  un 
large  intervalle  de  longueurs  d’onde  connues  donnent  des  informa¬ 
tions  ddfinies  sur  le  rayon  de  I’orbite  electronique  4f.  Nos  mesu 
res  indiquant  un  rayon  different  de  celui  employe  par  Blume  , 
Freeman  et  Watson,  qui  ont  calcuie  le  facteur  de  forme  en  assumant 
des  fonctions  d‘ onde  hydrogeniques .  II  sera  discute  le  choix 
possible  de  differentes  constantes  de  screening. 

ABSTRACT 

The  total  neutron  cross  section  for  holmium  was  measur 
ed  for  neutron  wavelengths  between  0*20  and  9«00  Angstroms.  The 
IEA  crystal  spectrometer  was  used  together  with  a  mechanical 
velocity  selector  to  eliminate  higher  order  reflections  from 
the  crystal.  Powder  samples  of  HogO-j  were  prepared  in  high 
purity,  talcing  particular  care  to  eliminate  contamination  by  the 
rare  earths  of  high  neutron  cross  section.  Although  paramagnetic 
form  factors  are  usually  determined  by  differential  cross  sector 
measurements  with  one  neutron  wave length,  aocurate  total  cross 
section  measurements  for  a  wide  range  of  known  wavelengths  give 
definite  information  on  the  radius  of  the  4f  electron  orbits  . 
Our  measurements  indicate  a  radius  different  from  the  one  used 
by  Blume,  Freeman  and  Watson,  who  calculated  the  form  factor 
assuming  hydrogenic  wave  functions.  The  possible  choise  of 
different  screening  constant  will  be  discussed. 
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I.  INTRODUCTION 


During  a  program  of  total  cross  section  measurements  of 
the  rare  earth  elements  at  this  Institute*  holmium  was  given 
special  attention  because  of  the  interesting  interaction  between 
the  neutrons  and  the  atomic  electrons.  This  paramagnetic 
scattering  has  been  noticed  before  in  total  neutron  cross 
seotions  by  Bernstein  jet  al .  ~  and  in  neutron  diffraction 
work  by  Koehler*  Wollan  and  Wilkinson2. 


The  total  neutron  cross  section  for  scattering  of  neutrons 
by  an  atom  in  the  presence  only  of  an  interaction  between  the 
magnetic  moment  of  the  neutron*  y  *  and  that  of  the  atom,^  »  is 
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where  — is  the  classical  electron'  radius,  and  f  is  the 
me*1 

neutron  scattering  form  factor  for  those  unpaired  electrons 
whioh  contribute  to  the  magnetic  moment  of  the  holmium  ion. 


The  magnetic  moment  of  the  holmium  ion  has  been  measured  as 
being  10.^4  1  olO  Bohr  magnetons  by  Strandburg*  Legvold  and 
Spedding^*  The  paramagnetic  scattering  fora  factor  has  been 
calculated  by  Blume,  Freeman  and  Watson^*  using  Hartree  ~  Fock 
wave  functions  for  isolated  ions.  These  form  factors  ,  together 
with  known  values  of  the  neutron  and  ionic  magnetic  moments,  may 
be  used  to  predict  the  paramagnetic  scattering  component  to  the 
neutron  scattering. 


HowevC"".  it  is  essential  to  make  accurate  measurements  of 
the  neutron  cross  section  to  supply  information  unavailable  from 
the  theory.  In  particular*  the  ionic  fora  factor  for  the  ion  in 
the  solid  state  may  be  expected  be  different  from  the  one  for 
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a  free  ion.  In  addition,  measurements  over  a  wide  range  of 
neutron  energies  are  required  to  determine  the  variations  of  the 
nuclear  cross  sections  with  neutron  energy. 

Ordinarily,  scattering  form  factors  are  best  determined 
<  by  differential  scattering  experiments.  Holmium  has  been  measur¬ 

ed  in  this  way  by  Koehler,  Wollan  and  Wilkinson2.  In  the  present 
»  work,  it  is  shown  that  careful  total  oross  section  measurements, 

besides  giving  the  nuclear  oross  sections,  give  information  about 
the  4f  electron  orbits. 

II.  EXPERIMENTAL  METHOD 

The  source  of  neutrons  for  this  work  was  the  Instituto  de 
Energia  A tunica  swimming  pool  research  reactor  operated  at  2  Mw. 
The  thermal  neutron  flux  measured  with  gold  foils  near  the  core 
of  the  reactor  is  2.4  x  1012  neutrons/cm2/seo. 

A  crystal  spectrometer  and  a  mechanical  velocity  selector, 
constructed  in  the  shops  of  the  Instituto  de  Energia  Atdtaioa  , 
were  used  as  monochromators .  Order  contamination  was  eliminated 
from  the  beam  by  using  the  orystal  together  with  the  mechanical 
velocity  seleotor. 

The  crystal  spectrometer  was  first  located  olose  to  the 
,  reactor  in  a  radial  beam-hole,  as  shown  in  Figure  1.  Later,  the 

crystal  and  the  selector  were  looated  at  one  of  the  tangential 
#  beam- holes,  as  shown  in  Figure  2.  We  had  a  thermal  flux  of 

4  x  10^  neutrons/cm2/seo.  outside  the  radial  port  and  of 

1.4  x  10^  neutrons/cm2/seo.  outside  the  tangential  port.  The 

crystal  spectrometer  angles  were  read  on  a  vernier  scale  with  a 
precision  of  0.01  degree.  The  measurements  were  made  with 
crystals  of  calcite,  aluminium  and  mica.  With  the  calcite 
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crystal  the  resolution  was  18  minutes. 

The  mechanical  velocity  selector5  gives  a  neutron  energy 
with  a  resolution  which  depends  on  the  rotation  velooity,  end  on 
the  helical  channel  inclination#  We  used  a  resolution  of  50^ 
in  wavelength#  The  final  resolution  was  determined  by  the 
crystal  spectrometer  and  order  contamination  was  complete¬ 
ly  eliminated  by  the  mechanical  velocity  selector# 

Commercial  boron  tri-fluoride  detector*,  enriched  in  the 
isotope  B10,  were  used  fee*  neutron  detection# 


Mterwt 
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Figure  1  »  Soheefttle  diagram  of  the  crystal 
spectre©®  ter 

The  samples  of  holmium  eo&da  were  supplied  by  the  Chemical 
Bagineering  Division  of  this  institute.  The  separation  method 
employed,  using  ion  exchange  resins,  assured  u»  of  the  degree 
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6* 

of  purification  nsoeasary  far  this  exporisttnte 

The  sssples  were  pi «w;,  in  iluadaim  containers  and 

introduced  Into  the  bean  in  a  reproducible  position®  Baa 

transjoiBBian  through  the  ssogple  was  obtained  by  »»aauring  the 

counting  rate  with  the  Maple  in  the  bean*  and  the  rate  obtained 

with  an  identloal  e^arty  ssc$le  holder  in  the  bean«  A  background 

was  subtracted  frora  each  uo«iiiag.  The  containers  were  designed 

to  give  a  tr&osaiision  which  idaisdaed  the  tias  required  to 

c 

reduce  the  statistical  errors  ® 


Figure  2  -  Schematic  diagreai  of  ths  mohanical 
velocity  selector  together  with  the  crystal 
speotroastsr* 

To  avoid  ths  influence  on  the  trsasaisaion  of  the  fluctua¬ 
tions  of  the  reactor  power  and  ths  instability  of  the 
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electronic  circuits,  the  traasmiasim  ©s&suressents  were  repeated 
several  tdma  in  cycles,  &ooord,ing  to  a  routine  designed  to  cancel 
1  inoar  drifts*  detector  puls©a  *?sre  ®ggpXifi©d,  analysed  and 

oomted  by  two  independent  elsotroale  system . 

total  oroos  sections  mr®  calculated  fro®  the 
transmission  B5eamtr«sant»*  2ba  oceventicml  formula  haa  bean 
used  for  the  error*  Th®  calculations  were  jsads  by  on  X£@f>l620 
oojaputer.  The  correction  due  to  oxyg«n  was  mde  siseply  by 
subtracting  the  free  atoa  ocxygen  cross  section  of  3*8  barn®  per 
a  tow.  Figure  3  shows  the  result*  in  the  usual  way  with  total 
cross  section  versus  heutron  energy* 


Figure  4  -  Total  neutron  oross  section  of  holmiuaj 
«®  a  function  of  neutron  wavelength.  The  oaloulat 
ed  paremgaatlo  areas  section,  ^pm  *  is  shown.  ~* 

m.  as&lxsxs  of  m 

'Brio  total  cross  section  of  holdu®,  (f^  ,  oonaiats  of  three 

completely  Independent  partial  crocs  sections,  one  of  which  Is 

the  par&asegaetlo  scattering,  (T  .  Than, 

pm 

. 
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wnere  O'*  ^  fl£  are  the  nuclear  scattering  and  absorption  . 
respectively.  was  assumed  to  be  independent  of  energy  and 
was  assumed  to  vary  as  1/^E  over  the  range  of  this  experiment. 
The  nuclear  resonance  spacing  is  about  6  electron  volts?.  The 
first  resonance  at  3»92  electron  volts  has  no  effect  on  these 
assumptions.  It  is  unlikely  that  a  bound  state  would  influence 
the  dependence  on  energy  of  the  nuclear  cross  sections. 

Nuclear  scattering  contributes  least  to  our  experimental 
data.  For  the  purpose  of  analysis*  was  taken  as  7  barns  by 
comparing  with  neighbouring  nuclei?.  Nuclear  absorption 
dominates  at  very  low  energies  and  may  be  quite  well  determined 
by  using  the  asymptotic  value  of  (T^  =  65.2  -  1.3  barns*  caloxilflt 
ed  from  the  holmiura  ion  magnetic  moment. 

The  nuclear  absorption  cross  section  thus  determined  from 
our  data  is  ^  *=  6l  -  3  bams,  reduoed  to  its  value  at  thermal 
neutron  energy  .025  eleotron  volts. 

IV.  CONCLUSIONS 

Figure  5  shows  our  experimental  points  after  having 

subtracted  the  contributions  of  nuclear  scattering  and  absorption. 

Comparing  our  experimentally  determined  paramagnetic 
scattering  cross  section  with  that  expected  from  the  Hartree- 
-Fock  calculations  of  Blume,  Freeman  and  Watson ^  one  sees  that 
the  disagreement  is  not  serious,  A  4f  shell  radius  smaller  by 
10 J3,  or  a  screening  constant  5#  less,  would  minimize  the 
discrepancy  between  their  theory  and  our  results.  However,  an 
agreement  within  experimental  errors  cannot  be  obtained  only  by 
changing  the  radius,  but  rather  by  admitting  the  possibility  of 
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a  slightly  mors  diffuse  wavo  function*  such  as  pro£x>s»d  toy  Judd 
and  Liadgron®.  Unfortunately*  no  paramagnetic  neutron  scattar¬ 
ing  fora  factors  are  available  using  thair  K*ve  functions. 


Figure  5  -  Comparison  between  the  axperirontal 
paramagnetic  cross  section*  obtained  after 
having  subtracted  the  contribution  of  nuclear 
scattering  and  absorption*  wad  the  theoretioal 
parasagnetlo  crocs  section*  calculated  from  the 
fore  factors  of  Slusse*  Fresoan  and  Watson4.  the 
older  theory  quoted  by  Koehler  _#t  &l$  is  shown  . 

The  apparently  discrepant  points  due  to  Beenatain 
et  alt  are  explained  In  the  text. 

Bie  determined  values  of  (J^  of  Bernstein  jst  al.  are 
different  from  ours  although  Figure  3  shows  that  the  total  cross 
sections  agree  In  the  region  where  their  statistical  error*  were 
comparable  with  ours.  Thair  paraaognetio  cross  sections  differs 
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from  our  values  which  have  taken  advantage  of  more  recent 
available  information.  Their  results  were  corrected  using  a 
nuclear  scattering  cross  seotion  of  13  barns,  based  on  a  total 
cross  section  of  28  bams  at  05  eV.  It  now  appears  from  the 
compilation  of  all  total  cross  section  values  that  28  bams  is 
at  least  3  bams  too  high,  as  it  can  be  seen  in  Figure  3*  In 
addition,  they  negleoted  the  effects  of  the  then  unknown  3*92  eV 
resonance  on  the  capture  cross  section  at  C£  eV  (about  1  bam), 
and  they  underestimated  (by  about  1  bam)  -die  paramagnetic  scatter¬ 
ing  cross  section  at  05  eV.  Taken  together,  their  estimate  of 
13  bams  for  nuclear  scattering  is  about  5  cams  too  high.  In 
our  lower  energy  range,  the  nuclear  scattering  is  relatively  un¬ 
important;  in  the  absence  of  direct  measurements  ,  we  used  7 
bams  based  on  measured  values  of  nuclear  radii  of  neighbouring 
nuclei1®*  Finally,  analysis  with  our  new  value  of  the  thermal 
absorption  cross  section  of  6l  -  3  bams  instead  of  64  assumed 
by  Bernstein  _et  al.  accounts  for  the  remaining  small  discrepancy 
between  the  analysis  of  our  results  and  that  of  their  careful 
measurements . 

Figure  5  also  shows  the  disagreement  between  the  older 

? 

calculation  quoted  in  the  experimental  work  of  Koehler  et  al., 

who  used  the  theory  by  Trammell11  with  hydrogen-like  wave  functions, 

a 

and  the  calculation  of  Blume,  Freeman  and  Watson  who  used  the 
same  theory  with  Hartree-Fock  wave  functions®  Blume,  Freeman  and 
Watson  noted  this  disagreement.  The  interpretation  of  our 
results,  according  to  Trammell's  theory,  favours  the  smaller  4f 
shell  implied  by  the  more  recent  calculations  of  Blume,  Freeman 


We  plan  to  continue  our  work  with  metallic  samples  to 
eliminate  small  uncertainties  in  subtracting  the  scattering  by 
oxigen.  In  addition,  it  is  planned  to  make  an  independent 
measurement  of  the  nuclear  absorption  cross  section  with 
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facilities  existing  in  this  Institute. 
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